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ABSTRACT 


This work deals with a study of the levels of the odd-odd nucleus 
64Cu, using the ®4Ni(p,ny)®*Cu reaction. Spin assignments, y-ray 
multipole mixing ratios and y-decay branching ratios have been 
determined for the levels below |] MeV excitation energy. The 
experimental level scheme is also compared with the predictions of the 
Intermediate Coupling model. 

Properties of the low-lying excited states of ®*Cu were 
determined by measuring the angular distributions, linear polarization 
angular distributions and angular triple correlations of the y rays 
emitted by these levels. Spin assignments and y-ray multipole 
mixing ratios were determined by comparing the experimental data with 
the predictions of the Statistical Compound Nuclear model. The 
phase of the mixing ratio employed in this study is consistent with 
that of Rose and Brink; an extension of their phase consistent 
formalism to include linear polarization angular distributions and 
y - y triple angular correlations is presented in an Appendix. 

The following spin-parity assignments have been deduced for 
levels in ©4Cu: 159 keV (J = 2), 278 keV (U" = o(t)) 344 keV 


(J” = ] Cl 


), 362 keV (J = 3), 574 keV (J = 4(2)), 609 kev (J" = 2 


> 


GOcekeVe eam limes OekeVa main? 3))) sme hoekeVe( duis) Mog oekevini d= —0) 


895 Kev (J 


| 


3) and 927 keV (J = 1(3)).  Multipole mixing ratios have 
also been determined for most of the y decays between these levels. 
The results obtained in this experiment are compared with previous 


works. 
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The properties of the low-lying energy levels of °*Cu have been 
calculated using the Intermediate Coupling model. The results of 
these calculations are in reasonable agreement with the excitation 
energies and single particle spectroscopic factors for many of the 
levels below 1 MeV excitation, although there are some aspects of 
the experimental spectrum which are not well reproduced. Suggestions 


are made for possible improvements to the calculation. 
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11. x? curves for the 278 > 0 angular distribution, 
showing the dependence of the angular 
distribution on the parity of the initial 
State in the decay. 


12. The 159 keV level: (a) x? curves and (b) angular 
distribution for the 159 > 0 transition. 

The 278 keV level: (c) x* curves and (d) angular 
distribution for the 270,~> 0, transitJwon. 


13. The 278 keV level: (a) y? curves and (b) linear 
polarization distribution for the 278 keV y ray. 
The solid lines in (b) show the expected polariz- 
ation distributions for values of J" and § which 
gave good fits to the angular distribution 


14. The 344 keV level: (a) x? curves, (b) angular 
distribution for the 185 keV y ray, (c) Nap 
geometry and (d) Cl geometry of the 344 + 159 + 0 
angular correlation. 


15. The 344 keV level: (a) x? curves and (b) angular 
distribution for the 344 + 0 transition. The 
362 keV level: (c) y? curves and (d) angular 
distribution for the 362 + 159 transition. 


16. The 344 keV level: (a) x? curves and (b) linear 
polarization distribution for the 344 keV 
y ray. 


17. The 574 keV level: (a) x? curves, (b) angular 
distribution for the 212 keV y ray, (c) Al 
geometry and (d) Cl geometry of the 574 > 362 > 159 
angular correlation. 


18. The 609 keV level: (a) x? curves, (b) angular 
distribution for the 450 keV transition, (c) Al 
geometry and (d) Cl geometry of the 609 + 159 + 0 
angular correlation. 


19. The 609 keV level: (a) y? curves and (b) angular 
distribution for the 609 keV y ray. (c) x? curves 
and (d) angular distribution for the 609 > 344 
transition. For (c) and (d), J refers to the 
spin of the 344 keV level, 


20. The 609 KeV level: (a) y* curves and (b) linear 
polarization distribution for the 609 keV y ray. 
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21. The 663 keV level: (a) y* curves, (b) angular 
distribution for the 385 keV transition, (c) Al 
geometry and (d) Cl geometry for the 663 + 278 + 0 
angular correlation. 61 


22. The 663 keV level: (a) x* curves and (b) angular 


distribution for the 663 keV transition. Gey 
curves and (d) angular distribution for the 
663 > 344 transition. 62 


23. The 739 keV level: (a) x? curves, (b) angular 
distribution for the 580 keV y ray (c) Al geometry 
and (d) Cl geometry of the 739 > 159 > 0 angular 
correlation. 65 


24. The 739 keV level: (a) xX? curves, (b) angular 
distribution for the 461 keV y ray (c) Al geometry 
and (d) Cl geometry of the 739 > 278 > 0 angular 
correlation. 66 


25. The 739 keV level: (a) xX? curves and (b) angular 
distribution for the 739 > 362 transition. 67 


26. The 746 keV level: (a) x* curves, (b) angular 
distribution for the 468 keV y ray (c) Al geometry 
and (d) Cl geometry of the 746 > 278 > 0 angular 
correlation. 68 


27. The 878 keV level: (a) x? curves and (b) angular 
distribution for the 878 > 0 transition. (c) x’ 
curves and (d) angular distribution for the 
878 + 159 transition. 70 


28. The 878 keV level: (a) x* curves, (b) angular 
distribution for the 534 keV Y ray, (c) Al 
geometry and (d) Cl geometry of the 878 > 344 + 0 
angular correlation. 7) 


29. The 895 keV level: (a) X* curves and (b) angular 
distribution for the 895 > 278 transition. (c) x 
curves and (d) angular distribution for the 
895 > 159 transition. 72 


30. The 927 keV level: (a) X* curves, (b) angular 
distribution for the 649 keV Y ray (c) Al 
geometry and (d) Cl geometry of the 927 * 278 > 0 
angular correlation. 73 
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CHAPTER 1 
MOTIVATION 


The odd-odd nucleus °*Cu has been the subject of several 
experimental investigations in recent years. The (d,p) and (d,) 
measurements of Park and Daehnick (Pa 69) established the level 
scheme up to about 3 MeV excitation, and measurement of the charged 
particle angular distributions resulted in the assignment of af! 
limits for most of the low-lying excited states. The (n,y) capture 
work of Shera and Bolotin (Sh 68) Doren ened the y-decay scheme for 
the levels below 1322 keV excitation in °*Cu, and provided tentative 
spin assignments for most of these levels based on the assumption 
that the predominant mode of decay is via MI radiation. 

Several authors (Ve 61, Mi 62, Ko 65, Ko 69) have carried out 
measurements on the circular polarization of the y rays from °*Cu 
following the capture of polarized thermal neutrons by °©%Cu. These 
measurements generally did not lead to unique spin assignments, and 
often produced results which were in serious disagreement with other 
works. 

In the last few years, several authors have exploited the 
S4Ni(p,ny)®*Cu reaction to study the low-lying levels of °*Cu. 
Davidson et al (Da 70) and Drum (Dr 70) have studied the double 
correlation (often called the angular distribution) of the y rays from 


levels of ©*Cu populated in the °*Ni(p,n) reaction, and have presented 


' 7 : - : - 
; = a: a A¢ 


- - Vv _ 7 i eT. ; _ 
; hg SE SY 
| ; Wauaeet if a ee 


letsyae To ide ia ont aad es 0" 7 aues ane 

{®,b) bane (oyb) it gay 3719997 At Seeotbaisaeiht fesmeaivegge 

+ (avel sda badgifdeses (03 6%) dstarsst® we wet to einemeiwessia 
bepists alt to tHsmenvesom brs \naissnidxs Ver £ su0ds oF qu smedoe 
"L Yo snamapizes adj ni badfu2er enolgudiizetb velupns slatsieq 
siuiqe2 (¢,n) ofT .2e%6te betioxs pniyi-wet ais to seem 70? axial! 


4ot émodoe yavab-y $14 badeffdstes® (B83 dz) qisolod bas eredé To dow a 
= ¢ ; } 


avivetus2 babivorq brs ,vd"? wi noltezioxs Ved SSEf wolsd efovel ods 
hoizamuers ody no bezed efavel szeris to t20m sot 2inemnpizes niqe : 

. sndidstbsy (M siv 2l ysosb Ye sbom tneaimoabsig sft tenz 

tuo beinisa syed (23 of (28 oF ;83 1M ,13 9V) 210nsuB lereve? 
ud"* mort 2yet y Sit Io nolsssiialoq tefustin sd? po stnemetueeem 
eesti =. F?_s vd ancysusn (swiadt -bositeiog Yo swwideo edz pniwoffo? 7 
bre ,2snemnpieee nige supinu o¢ bas! ton bib yffersnep etnsmewesom | 
jafto diiw inemsstoselb eupits2 ni siaw doifw 2tfuess beouborg neato 
.2a10W 

sii batliolqzs eved 210ornsus levevee 2isay ws? tee! eda al ai 

ud"? Io etavel oniyi-wol srf1 ybute of noitgees va? bye) HO 

- gldvob saz bothuse eved | (G§ 790) mpuvd dns (OT 5a) isis minhaa 

| mont ayer y sda Yo (nOTudirieTb yelupne edt bal fsa agate) vies 
basnezeng over bns enincade (n,q) tHe? =F ene we ihe oo 


spin assignments based on the comparison of the data with the 
predictions of the Statistical Compound Nuclear model (Sh 63, Sh 65). 
In several cases, it was not possible to make unique spin assignments 
on the basis of the angular distribution alone (see Section 2.2.1), 
and for many of the y decays, the multipole mixing ratio could not be 
uniquely determined. While the results of these two studies are 
generally consistent with each other, several spin assignments do not 
agree with previous works. Wel lborn et_al (We 71) have studied y - y 
triple correlations from ©*Cu, again using the ®*Ni(p,n) reaction to 
excite the states of interest. Their results are in agreement with 
those of Davidson et_al (Da 70) and Drum (Dr 70), but this work dealt 
mainly with the states below 400 keV excitation. Mischenko et al 
(Mi 72) have studied (n,y) correlations following the ©*Ni(p,ny)®*Cu 
reaction, and have deduced spin assignments for the first four excited 
States Of .° GU. 

In spite of the wealth of experimental data available for the 
low-lying states of ®*Cu, the situation regarding spin assignments 
for these levels, especially the levels above 400 keV excitation, and 
multipole mixing ratios for the y decays is still somewhat uncertain. 
For this reason, we decided to re-investigate the y decay of this 
nucleus, using a combination of double correlation, linear polarization 
double correlation and y - y triple correlation measurements to 
attempt to resolve the present ambiguities. 

We have also chosen the °*Ni(p,n) reaction to excite this nucleus, 
for several reasons. The relatively large cross-section for this 


reaction will permit measurements to be performed for many of the 
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weaker y transitions, and will] provide a high enough yield for the 
coincidence experiments to be feasible. In this region of the 
periodic table, the mechanism of the (p,n) reaction can be described 
by the Statistical Compound Nuclear Theory, which can be used to 
calculate the substate population parameters of the y-decaying level 
provided that this state is not populated by a cascade transition from 
a higher excited state. The endothermic nature of the reaction 

(Q = -2.461 MeV) allows the state of interest to be excited close to 
threshold, thereby eliminating complications due to cascade feeding. 
Moreover, the usual target contaminants, carbon and oxygen, have very 
large negative Q-values for protons, which results in a fairly ''clean'' 
y-ray Spectrum. 

From a theoretical standpoint, there have been no attempts to 
apply any of the usual nuclear models to ®*Cu. Shell model 
calculations, in particular, would be exceedingly difficult, due to 
the prohibitively large configuration space which would be necessary. 
On the other hand, there is no evidence that nuclei in this region of 
the periodic table are permanently deformed, so Nilsson model 
calculations have little chance of success. 

The most promising approach to a theoretical study of ©*Cu is 
probably via the Intermediate Coupling model (Ch 54, Ch 67). The 
even-even nucleus °*Ni shows a characteristic vibrational-type 
spectrum, and both Gomez (Go 71) and Markham (Ma 71) have given a 


description of ©3Cu as a single proton coupled to this vibrating core. 
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The success of these calculations for °*Cu suggests that it might be 
possible to describe low-lying states in °*Cu in terms of two particles 


(a proton and a neutron) coupled to the surface vibrations of °?Ni. 
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CHAPTER 2 


THEORY 


2.) Phase Consistency in_y Decay 


The technique of observing the angular distribution or angular 
triple correlation of the y rays emitted by excited nuclear states is 
a widely used source of information on the properties of these states. 
In particular, such measurements can be used to determine the 
multipole mixing ratios (i.e. ratios of reduced matrix elements for 
different multipoles) of the decay, which are sensitive to the 
detailed structure of the nuclear wave functions. 

The theory of y-ray angular correlations has been presented by 
many authors (Bi 53, De 57, Bi 60, Fe 65, Po 65), but these earlier 
formulations all suffered from the lack of an explicit definition of 
the phase of the mixing ratio. This in turn made it extremely 
difficult to compare the sign of a measured mixing ratio with the 
predictions of a nuclear theory. This problem has been overcome with 
the development of the theory by Rose and Brink (Ro 67), who adopted 
the convention that the electromagnetic transition multipole operators 


<1> ; ; 
Tis T should transform under time reversal according to the rule 
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where 6 is the time reversal operator. 

The work of Rose and Brink (Ro 67) included only the development 
of the angular distribution formula. [In order to make use of their 
phase consistent theory in the present work, their formalism has been 
extended to include linear polarization and angular triple correlations 
as well (Gr 72). This development is given in Appendix A. 

The next section contains a brief discussion of the three 
different measurements which were performed in this experiment, with 
particular emphasis on what new information can be obtained from each. 


A detailed derivation of the formulae used can be found in Appendix A. 


2.2 Angular Correlation Theory 


2.2.1 y-Ray Angular Distribution (Double Correlation) 
The angular distribution of a y ray refers to the relative 
intensity of the y ray yield as a function of angle, and can be 


expressed in the form 


cos 6) (2-2) 


W(0) = Py 
geakaek 


where 6 is the angle between the direction of emission and the 
z-axis. The a, coefficients in (2-2) depend on the spins of the 
initial and final nuclear states, the relative amplitudes of the 
different multipoles contributing to the decay, and the mechanism of 
the reaction leading to the formation of the initial state. 

The amount of useful information which can be extracted from an 


angular distribution measurement is governed by the number of ay 


coefficients present in the expansion (2-2), and this is limited by 
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the following considerations. 
i) If the initial and final states are states of 

definite parity, and if parity is conserved by the 

electromagnetic interaction, then the summation over 

k is limited to even values only. 

ii) The subscript k must satisfy the triangle 

conditions A(L',L',k) and A(J,,J,,k), where L' is the 

highest multipole in the transition, and J, is the spin 

of the initial state. |The symbol A(a,b,c) means 

faa p< eCeceat. 

Except in very rare cases (see e.g. Li 64), it is a good 
assumption than multipoles with L > 3 do not contribute significantly 
to the decay. It is therefore possible to extract, at most, the 
a, and a, coefficient from the data (a) is just a normalization 
factor). Usually, the final spin in the decay is known, so the 
unknown parameters are the initial spin, a single quadrupole/dipole 
mixing ratio (6), and the substate population parameters. 

The angular distribution formula (2-2) assumes that the radiation 
is observed by a point detector. In practice, detectors always have 
a finite angular spread, and the result of averaging the angular 
distribution over the solid angle of the detector is to reduce the 
measured anisotropy of the distribution. Provided that the detector 
has cylindrical symmetry, one can correct for this effect exactly by 
multiplying each term in (2-2) by an "'attenuation coefficient'' 


Q (Ro 53), and the measured angular distribution is 
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W. (6) = : a, Q P, (cos 9) (2-3) 


In certain instances, it is possible to adjust the conditions of 
the experiment so that the population parameters can be calculated 
from a reaction theory. In the present work, the Statistical Compound 
Nuclear model has been used, and a brief discussion of the assumptions 
of this model will be presented in Section 2.4. Even with this 
additional tool, however, it is often impossible to uniquely determine 
both the spin and mixing ratio from an angular distribution measurement 
alone, and it is then necessary to resort to other techniques, such as 
linear polarization or triple correlation measurements in order to 


resolve these ambiguities. 


2.2.2 y-Ray Linear Polarization Double Correlation 
The electromagnetic multipole fields of order L can be divided 
into two classes, the electric multipoles (with parity Ge and the 


ake At large distances from 


magnetic multipoles (with parity (- 
the source of the radiation, the only difference between the electric 
and magnetic multipole fields is that the electric and magnetic vectors 
of the radiation are both rotated through an angle of 1/2 about the 
direction of propagation. In order to measure the parity of the 
radiation, it is therefore necessary to determine the direction of 
the electric (or magnetic) vector. 

The linear polarization of a y ray is defined as the direction of 


the electric vector. If @ is the angle between the reaction plane 


(the plane containing both the z-axis and the y ray) and the 
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polarization direction of the radiation, the angular distribution 


of linearly polarized photons is given by 


I easrsvony ld uaa ait Vea 


W(6,¢) = r{a, P, (cos ovat (-) a 
k 


where the index 7' is zero if the L' radiation is electric and one if 


it is magnetic. The degree of linear polarization is defined as 


since this limits p(6) to lie between +1 and -1, with p(6) = 0 for an 
unpolarized y ray. 

Apart from the dependence of the linear polarization double 
correlation on the parity of the radiation (and hence on the relative 
parities of the initial and final nuclear states), the measurement of 


the coefficients a' in (2-5) can give further information on the spins 


k 
and mixing ratio. For example, there are certain values of 6 for 
which the angular distribution will be isotropic, but the y ray may 


nevertheless be strongly polarized (Gr 71), and a measurement of p(6) 


can eliminate some possible spin assignments. 


2.2.3 y-Ray Triple Correlations 
A third method of obtaining additional information about the 
decay is through the measurement of two y rays in coincidence. The 


vy - y triple correlation can be expressed in the form 
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W(6, »O ») = Me (2-6) 
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and, as before, the a coefficients are functions of the spins and 


ki kok 
multipoles involved in each of the decays, and also of the population 
parameters of the initial state. The subscripts k,, kp and k3 are 
limited to even values only, and are also governed by the restrictions 
ACS sky) secily glaun ks) yo AdLetslo shad A(do,do,ks)eand Aik: ,koyks) 
If multipoles higher than quadrupole do not contribute to either decay, 
there are a maximum of 19 coefficients which can be determined 
ECan Mental IVA and although the number of unknown parameters is 
increased, there are enough measurements possible to overdetermine the 
problem. 

In practice, the number of unknowns in the problem is unlikely 
to exceed about five, and it is therefore unnecessary to measure all of 
the 19 coefficients in equation (2-6). Instead, it is usually 
sufficient to determine certain linear combinations of these 
coefficients, by measuring the correlation in one or more ''geometries", 
in which two of the three angles 91, 92, > are held fixed, and the 
third is allowed to vary. lf the fixed angles are chosen properly, 
it is possible to write the correlation as an ordinary Legendre 


polynomial series, by using the expansion 


k 
6,,0 = 2 S) - 
cakikes 1,02,9) k Bkikoks Py (cos ) (2-7) 
where 9 is the single variable angle. The index k takes on even values 


only, and is limited to 0 < k < 4 if multipoles higher than quadrupole 
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are not included. Finite geometry effects are taken into account by 
multiplying each term in equation (2-6) by the factor Q. Q, , where 
2 Kg 
Q. is the attenuation coefficient for the detector of the first y ray 
2 


in the cascade, and Q. is the attenuation coefficient for the detector 
3 


of the second radiation. 


2.3 Measurement of Linear Polarization Using Compton Scattering 


In order to measure the linear polarization of ay ray it is 
necessary to make use of a process which is sensitive to the direction 
of the electric vector of the radiation. Although several processes 
have been suggested (Fa 59), the most commonly used has been Compton 
Scattering, because of the relatively large cross-section and the high 
sensitivity to linear polarization over an energy range containing 
most nuclear y rays of interest. 

In practice, the polarization of the scattered quantum is not 
observed, and the differential cross-section for the Compton scattering 


of a linearly polarized photon is (Fa 59) 


dopmt penn (Ky * (Koy Ka oo sin? 2 2 
ae ie is oe 2esin-Wicos7 7c) (2-8) 


where r, = e*/m,c? is the classical radius of the electron, yp is the 
Compton scattering angle, and xz is the angle between the polarization 
direction of the incident radiation and the scattering plane. ko is 
the energy (in units of m)c*) of the incident radiation, and 


k = kp/{] + ky (1 - cos w)} is the energy of the scattered photon. 
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The linear polarization p(6) is measured by observing the number 
of quanta Ng and Ngg which are scattered parallel to the reaction plane 
and perpendicular to the reaction plane respectively. Neglecting 


solid angle effects, these numbers are given by (Tw 70) 


Ny = W(0,0) dog + W(8,90) dog, 


Noo W(6,0) dogo + W(9,90) do 9 (2-9) 


where ddg and dOgp are respectively the Compton cross-sections for 
scattering parallel to and perpendicular to the plane defined by the 
propagation direction of the incident y ray (ko) and its electric 
vector. The linear polarization (equation 2-5) can then be calculated 


from the expression 


p(e) = & Ge) 


Noo + No 
= d099 = dO 9 * sin*w a 
idoayHerdoy Ke Reis a {= 40) 
To 1 Sun wW 


where R, the asymmetry ratio, gives the response of the process to 
linear polarization. 

At low energies, the asymmetry ratio approaches unity for W = 90°, 
indicating that Compton scattering exhibits ideal response to linear 
polarization. Although R decreases quite rapidly with increasing 
Y-ray energy, Compton scattering has been used successfully for 
polarization measurements on Y rays of energy as high as 8 MeV (Su 59). 

As the photon energy increases, the ''best'' Compton scattering 
angle w decreases. Suffert et al (Su 59) have shown that, for 


MAX 


experiments which are limited primarily by counting statistics, the 
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best scattering angle is that which maximizes the quantity 


S (yp) Sek coe (2-11) 


where oF is the Compton cross-section averaged over polarization. 
Figure 1 shows a plot of this angle bax as a function of energy, as 
well as the asymmetry ratio R for yw = Day * 

The value of R given by equation (2-10) is an upper limit. In 
practice, the finite solid angles of the detectors used to measure the 


linear polarization will reduce the actual asymmetry R of a polarimeter. 


These effects will be discussed further in Section 3.3.1]. 


2.4 The Statistical Compound Nuclear Model 


The Statistical Compound Nuclear (SCN) model has been used with 
great success to describe the (p,ny) reaction on medium weight nuclei 
(Pi 70, Tw 70, Da 70, Mc 70,.Ro 71). Since the model has been treated 
in detail elsewhere (Sh 63, Sh 65, Vo 68), only a general outline will 
be given here. 

On quite general grounds, the differential cross-section for 


going from a reaction channel a to another channel a' is given by 


(Vo 68) 
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Figure 1. The Compton scattering angle Vy for which the 
asymmetry ratio R is a maximum, and the value of 


R for p=W,,,, plotted as a function of the energy 
MAX 
of the incident y ray. 
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where | is the target spin, i the projectile spin. Ss = |! + i| is 
the initial state channel spin, 2 is the relative orbital angular 
momentum of the two particles @ in the initial state and J = [s = a 
is the total angular momentum of the system. The primed variables 
refer to similar quantities in the exit channel, and the sum is over 
all of these indices, as well as the projections m— and mou of the 


channel spins. v2 


1 Many ees is an element of the collision matrix, 
> 


which connects the amplitude of the outgoing wave in channel a! to 
the incoming wave in channel a. The Z-coefficients contain only 
angular momentum coupling coefficients. 

-For compound nuclear processes, it is assumed that the (complex) 
numbers U vary rapidly with energy about a mean value of zero, and if 
we therefore average the cross-section over a sufficiently large 
energy interval, all terms linear in U will vanish. It is further 
assumed that the energy variations of different components of the 
scattering matrix are uncorrelated, and the only terms in (2-12) which 
survive are squares of the matrix elements. This is the statistical 
assumption. 

The compound nucleus assumption states that the decay of the 
compound nucleus is independent of its mode of formation, and hence 
each absolute square of a collision matrix element can be written as 
a formation probability multiplied by a decay branching ratio. 
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where T, _ ya) = 1 - je Sag}? is the optical model transmission 
coefficient. 

Consider specifically the (p,ny) reaction near threshold. Since 
the quantization axis is chosen along the beam direction, the projection 
of the incoming orbital angular momentum 2 is zero, and for the state 
J in the compound nucleus, only substates with M < I + 1/2 can be 
populated. Since the emerging neutron can be emitted in all 
directions, the projection of its orbital angular momentum is not 
necessarily zero, but if the neutrons are of low energy, &' = 0 partial 
waves will predominate, and only the contributions due to the particle 
spin can de-align the final state. For (p,ny) reactions near 
threshold, it is therefore expected that magnetic substates for 
M< | + 1 will be very nearly equally populated, and the populations 
of higher substates will be substantially smaller. lf all possible 
substates in the final nucleus are equally populated, the resulting 
y-ray angular distribution will be isotropic, and it is expected 


therefore that states with final spin < | + 1 will have isotropic 


decays. 


2.5 Data Analysis 

Spin-parity assignments and multipole mixing ratios are determined 
by comparing the experimental data with the predictions of equations 
(2-2), (2-5) and (2-6). <A quantitative measure of the ''goodness of 


fit'' is obtained by calculating the quantity 


S=y (W. - Y.)? (2-14) 
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where W is the theoretically predicted result for a measurement 

at the angle o., Y. is the value actually meacuredh and a.” is the 
variance in Yo. The measurements Y. have the Poisson distribution, 
which approximates the normal distribution if the Y. are large enough 
(say > 30). If this approximation is justified, then the minimum 
value of S has the x* distribution with N-p degrees of freedom, where 
N is the number of data points, and p the number of parameters which 
were varied to minimize S. For a single angular distribution 
measurement, p iS equal to 2, since an ao coefficient and a mixing 
ratio 6 are fitted to the data. 

Spin and mixing ratio combinations are accepted as possible 
solutions if the minimum value of the statistic S falls below the 
0.1% confidence limit of the xy? distribution. |The meaning of this 
confidence limit is that, if the theory used is the correct one 
(i.e. if the assumed spin and mixing ratio are correct), there is only 
one chance in 1,000 that S would be larger than this value. 

For each possible spin and mixing ratio combination, the error 
in the mixing ratio is then calculated using the 10% confidence limit 


of the F distribution, as suggested by Cline and Lesser (C1 70). 
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CHAPTER 3 


EXPERIMENTAL PROCEDURE 


3.1 General 


In this experiment, targets of Ni metal (enriched to 98% in ®*Ni) 
were prepared by combining the powdered metal with a mixture of 
polyurethane and benzene, and smearing the resulting glue onto a piece 
of 0.005'' Ta. Targets prepared in this manner were about 2 mg/cm? 
thick, and the resulting energy loss of ~ 100 keV for the incident 
protons was assumed to provide sufficient energy averaging over. 
compound nuclear resonances to allow the Statistical Compound Nuclear 
model to be used. 

The targets were bombarded with protons from the University of 
Alberta CN van de Graaff accelerator, with the beam collimated to a 
diameter of 1.5 mm at the target position. Beam currents were 
typically about 300 na. 

After a preliminary study of the yield curve to determine the 
decay scheme of the low-lying levels of ®*Cu, four bombarding energies 
were chosen for further measurements. These energies were 2.75, 3.00, 
3.30 and 3.75 MeV, and were chosen to eliminate complications arising 
from cascade feeding of the levels of interest from higher excited 
states. Many authors (Pi 69, Da 70, Dr 70) have stated that it is 


advisable to excite the y-decaying state near threshold, and hence 
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ensure that the exit neutrons are essentially pure s-wave, in order to 
maintain the alignment of the compound nuclear state. We have found, 
however, that the decrease in alignment at higher bombarding energies 

(say 500 keV above threshold) is not critical, and the increase in 


y-ray yield favours this procedure. 


3.2 Angular Distribution Measurements 


For the runs at 2.75, 3.00 and 3.30 MeV, the y rays were detected 
in a 38 cc Ge(Li) detector*, which had a resolution of typically 3 keV 
(FWHM) for the 1.332 MeV y ray from ®°Co. Pulses from the detector 
were fed through an amplifier**, and then into an analog to digital 
Ponverteul which was interfaced to a Honeywell DDP-516 computer. For 
the runs at 3.75 MeV, a 48 cc fetectorsa was used, which had a 
resolution of 2.5 keV (FWHM) at 1.332 MeV. 

For the runs at 3.30 and 3.75 MeV, measurements were taken at 


angles of 0°, 31°, 55°, 70° and 90°, since these angles are suitable 


for rapid extraction of the ay and a, coefficients from the experimental 


data (Gra 72a)e For the runs at the two lower energies, measurements 


° 


were taken at 6 = 45° instead of 9 = 55°. All angles were repeated 
once. Since measurements made at the same angle are statistically 


independent, they were used as separate points in the x? analysis of 


the data, and were not averaged to give a single point. 


“as Nucleare Diodes ,ilnc.gamodela LGCC 5.7 + 2:3 
** ORTEC, model 450 

+ TMC, model 217A 

++ ORTEC, model 8101-0724 
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A typical y-ray spectrum, taken at a bombarding energy of 3.75 
MeV, is shown in Figure 2. y rays corresponding to transitions 
in ©*Cu are labelled with the energies of the initial and final states 
in the decay, while contaminant peaks are indicated by the nucleus or 
the reaction responsible for the transition. The major contaminant 
peaks are the 18!Ta_ lines at 136, 165 and 301 keV, the *°Na y ray 
at 439 keV, the 74’72Ge(n,n') lines at 596 and 692 keV, and the 
S4Ni(p,yy)®5Cu peaks at 770 and 1115 keV. 

For the runs at 2.75, 3.00 and 3.30 MeV, the 770 keV y ray from 
the J = 1/2 first excited state of ®°Cu was used as an internal 
monitor. At 3.75 MeV, this transition could not be used since it was 
obscured by the 927 > 159 transition in ©*Cu, so one of the previously 
measured isotropic decays from a low-lying state of ®*Cu was used for 


an internal monitor. 


3.3 Linear Polarization Measurements 


3.3.1 Description of the Polarimeter 

The polarimeter used in this experiment was based on the Compton 
scattering process. In its simplest form, such a polarimeter consists 
of a Compton scatterer, and an analyzer to detect the Compton scattered 
radiation, but in practice, it is usual to also measure the energy 
deposited in the scatterer and hence reduce background by recording 
coincidences between the two detectors.. A Ge(Li) detector was used 
as a scatterer in the polarimeter, and a 3'' x 3!' Nal(TI)* detector 
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Harshaw ''Integral line assembly'' type 12S 12/3 
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A typical y-ray singles spectrum, taken at a proton 


bombarding energy of 3.75 MeV and an angle of 


Figure 2. 


55°, with the 48 cc. detector. 


§= 


was used as an analyzing crystal.  5'' of lead were used to shield the 
analyzing crystal from primary target radiation. The mean Compton 
scattering angle from scatterer to analyzer was chosen as wt = 82° 

Measurements of the polarization are usually made at two angles, 
Oe=)Ganandegs=_ 905. In order to reduce data collection time and 
errors due to repositioning of the analyzing crystal, the polarimeter 
was constructed with two analyzing detectors, so that the 0° and 90° 
measurements could be made simultaneously. A schematic diagram of the 
polarimeter is shown in Figure 3. 

For an ideal polarimeter (i.e. point scatterer and point 
analyzers), the asymmetry ratio is given by equation (2-10). In 
practice, the finite solid angles of the detectors must be taken 
into account, and the ideal value of R is never realized. The 
actual asymmetry ratio R for the polarimeter was calculated by 
numerical integration over the scattering and analyzing crystals, 
following the method outlined by Taras and Matas (Ta 68). The 
results of these calculations are shown in Figure 4. For energies 
above 1.5 MeV, the effect of the finite sizes of the detectors is 
small, but as the incident photon energy decreases, the error incurred 
by neglecting finite geometry effects becomes progressively larger, and 
for energies around 300 keV (which is the range we are interested in) 
can become as large as 30%. 

The actual asymmetry of the polarimeter was determined at two 
energies by measuring the polarization of the ®°Co y rays. The 
results of these measurements are also shown in Figure 4. The 


measured values agree with the calculated asymmetry to within 102, 
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Figure 3. A schematic diagram of the Compton. polarimeter used 
for the linear polarization measurements. 
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Figure 4. A comparison between the asymmetry ratio Ri EA for 
a Compton polarimeter consisting of point dePettors, 
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which was the experimental error in the measurements. 

The effects of two consecutive Compton scatterings have been 
ignored in this calculation of R. This effect is expected to be 
most important at lower incident photon energies, since the Compton 
scattering cross-section decreases for increasing photon energy. 
The agreement (within experimental error) between the calculated results 
and measured values of R for the ®°Co y rays indicates that this effect 
is probably small for incident energies greater than 1] MeV, but it is 
difficult to estimate the contribution from this process at lower 
photon energies. For the linear polarization measurements in ©*Cu, 


the calculated values of R were assigned an error of 10%. 


323.2 Linear Polarization Electronics 


A sum-coincidence circuit (Co 56, Ho 58) was used for each 
Nal-Ge(Li) pair, in which the parameter actually recorded is the 
sum of the energies deposited in the scattering and analyzing crystals. 
Recording the total energy for each event, rather than the energy 
deposited in one or the other of the detectors, has important 
advantages in the increased resolution obtained, especially if the 
analyzing crystal subtends a large angle at the scatterer. Even 
though, for a valid coincidence event, the energy deposited in one of 
the crystals varies with the Compton scattering angle, the sum of the 
energies is of course a constant, and the resolution is then limited 
only by the inherent resolution of the detectors (in this case, the 


Nal crystal). 
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A schematic diagram of the sum-coincidence electronics used in 
this experiment is shown in Figure 5. A fast timing signal for the 
Ge(Li) is obtained by feeding the output through a timing filter 
amplifier (TFA) and a constant fraction of pulse height trigger 
(CFPHT). This timing signal is delayed, and used to stop a time to 
amplitude converter (TAC). Start pulses for the TAC are obtained 
by feeding the output from the Nal through a double delay line 
amplifier (DDL) and a timing single channel analyzer (TSCA), which 
was also used to place an energy window on the scattered y rays. The 
TAC output is then fed through an SCA, to provide a coincidence gate 
for the sum spectrum. The timing resolution from this curcuit was 
about 20 ns FWHM. 

The portion of the electronics shown inside the dashed line was 
duplicated for the 90° Nal detector. As well as the two sum spectra 
and associated time spectra, the Ge(Li) singles spectrum was also 


recorded for the angular distribution measurements. 
3.3.3 Experimental Measurements 


Linear polarization measurements were carried out at a proton 
bombarding energy of 3.3 MeV. Measurements were taken at angles of 
Ome 3 bee 5 sen. and 90. . 

A typical sum-coincidence spectrum is shown in Figure 6. The 
spectra shown in this figure have been corrected for the different 
efficiencies of the two analyzing crystals by normalizing the results 
at @ = 0°, where the y rays are unpolarized. The peaks at 278, 344, 


and 609 keV are due to transitions in ®*Cu, and these were the only 
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A block diagram of the electronics used for the 
linear polarization measurements. The portion 
inside the dashed line was duplicated for the 
6=90° analyzing crystal. 
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Figure 6. A typical sum-coincidence spectrum, taken at a 
bombarding energy of 3.30 MeV and an angle of 
8=90°. 
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three transitions for which polarization measurements were possible. 
The peaks at 439 and 511 keV are due to the *°Na(p,p'y) reaction and 
Sy annihilation respectively. The strong 159 keV y ray in ®*Cu 
does not appear clearly in the sum spectrum, because of the poor 
timing qualities of the Ge(Li) detector at such low energies (the 


Compton scattered electron has an energy of only 34 keV), so it was 


not possible to perform a polarization measurement for this transition. 


3.4 y - y Angular Correlation Measurements 


3.4.1 Choice of Geometries 

As was mentioned in Section 2.2.3, it is usually unnecessary to 
measure all 19 coefficients of the y - y triple correlation, and one 
therefore measures the correlation in one or more ''geometries'', in 
which two out of the three angles specifying the correlation are held 
fixed. A list of the angles specifying the seven ''standard'' 
geometries of Ferguson (Fe 65) is given in Table |. inthis 
experiment, the Al and Cl geometries of the correlation were measured 
for each y ray cascade. 

The Cl geometry is especially useful for resolving ambiguities 
in spin and mixing ratio from the angular distribution measurements in 
ee CUmeeeLnestardet nucleuse> Ns has va Oo” ground state spin, and hence 
nearly isotropic y-ray angular distributions are expected for states 
with J = 0 or 1 (Section 2.4). For J = 2 states, there will be two 
specific values of § for which the a» coefficient vanishes (regardless 
of the degree of alignment), and if the corresponding ay, coefficient 


is small, these angular distributions will appear isotropic as well. 
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Table |! Angles specifying the seven y - y angular correlation 
geometries of Ferguson (Fe 65). 


Geometry ce 05 p 

Al Variable 90° 180° 
A2 90° Variable 180° 
BI Variable 0° 0° 
B2 On Variable 0° 
Cl Variable 90° 90° 
C2 90° Variable 90° 
D slo 90° Variable 


It is therefore expected that, if the angular distribution is isotropic, 
possible fits will be found for J = 0, 1 or 2. This behavior is 
clearly evident in Figure 15 for the 344 keV transition. For the 
decay of a spherically symmetric initial state, the angular 

correlation of two succeeding y rays depends only on the relative 

angle between the two y rays. In the Cl geometry, this relative 

angle is constant at 90°, and hence the Cl geometry will appear 
isotropic for a weakly aligned state. This measurement can then 


serve to accept or eliminate the J = 2 assignment. 


3.4.2 Electronics 

A schematic diagram of the electronics circuit for the angular 
correlation measurements is shown in Figure 7. The timing circuits 
for both the Ge(Li) and the Nal are identical with those described in 


Section 3.3.2 for the linear polarization measurements. The timing 
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resolution was again about 20 ns FWHM. 

The TAC output was fed through two SCA's, one corresponding 
to real + random coincidences, the other to random ceding tenes only. 
The SCA outputs were then used to generate routing signals for the 
Nal spectrum in ADC A, and were OR'd to give a coincidence gate for 
the Nal. The two Nal coincidence gates were then OR'd to give a 
general coincidence gate for ADC A and ADC B. 

Due to the organization of the computer software, it was 
necessary to put the signals from both Nal detectors into ADC A. The 
Nal analog signals were therefore first fed through a summing 
amplifier, which acted as an impedance matching device. The 
coincidence pulse for each Nal detector was used as an anticoincidence 
pulse to prevent the accidental summing of pulses from the two 
detectors. The portion of the circuit shown inside the dashed lines 


was repeated for the Nal detector in the Al geometry. 


3.4.3 Data Collection 

The on-line computer allowed sorting of the Ge(Li) spectrum in 
ADC B according to the ADC A spectrum. Four windows were set on 
each of the four Nal spectra in ADC A (one real + random and one 
random spectrum for each detector), giving sixteen spectra of y rays 
in coincidence with the 159, 278, 344 and 203 keV transitions from 
the first four excited states of ®'Cu. Runs were taken at Ge(Li) 
angles of 0°, 31°, 55°, 70° and 90°, all angles being repeated once. 


Each run was typically 5 - 6 hours long. 


— ytine gesrisbl anids mobiis? of 49 


ee at Se "wwe Sawn: 
b | 7 rt x 
| ae Mp uae 


unibnogear on bal +180? 


a oP 


as evant onto amen + leat 7 
ady 16? efenple onisuon aqeneneg'ay baeu merit sew etuqiuo ane ot 
not 3260 sonshisniea 6. svip ot b'A0 sxaw bre ,A 3GA nl minsosqge ten 
& SVip os b'AG nerd siew 2546p ssnsbionton let owt ofT 8.1 eH odd 
.@ 90A bhe A 20A 197 otee sarsbisnios fsvenag. : 
26w 1} ystewito2 Tesuqmos sit To foitesinspio eft o2 sud a 
gait .A IGA ofni esotoateb ish diod moat efengie sid Juq oF vragesoon 
pninmue ¢ Apuotit bet serit siotsisit aisw elenaizg polsns ish 
aft .saiveb enidozem sonsbeqni 5 26 betls datedw eititqns 
sonetiealoatine as 2o beeu 2aw ioyasteb bef dase 10? eetuq sonabi antes: 
ows sta mor? za2elug to enimmvé Istnebioos ond Ineveig 0 sefug 
vent! badesh odd sbieni nwore iludiis sfz To nolji0g sfflT . 2701 2838b 


.yiremosp fA ¢f3 ni rozastab lel sanz sot bsisesqs) 26w 


noitastfod e760 £.AcE . 
ni mintsaq2 (tJ)ed sat to onisio2 bawolls ssiuqmoz sait-no dT 
Ao Jee siaw ewoboiw qu04 .murtoege A 20A efy of gnibiosss @ 2A 
ano bre. fobaisy + leat so) A IGA ni svtasgqe Is wot oft to dose 
if. 


2y61 Y¥ to stizege nesdxie pnivip ,(ros2939b d5es FO? muys20qg2 mobs? | 


nor? enoitizenais Ved GOS bas ade .BYS , C21 ors da iw sonebTomios at 


(iJ)98 xe°naxtsd syow envA v0" 2o 2etst2 bedioxe wot sani? ban 


#90 baisegey pried esfons ile ,°0@ bas “OY .°22 .°1E .°O to zofer 
7 Pri 
-Rnol eyuon 3 - 2 ylisoiqy? caw wn 7 


33 


A typical coincidence spectrum is shown in Figure 8, showing 
y rays feeding the 278 keV level. This spectrum shows the Al geometry 
of the angular correlation, taken at a proton bombarding energy of 
3.75 MeV and an angle of 55°. The real coincidences with the 331] keV, 


385 keV, 461 keV, 468 keV and 649 keV transitions can be clearly seen. 


3.5 Data Analysis 


y-ray intensities were obtained directly from the computer display, 
using the program SUMMIT (Gr 71A), which fits either a linear or 
quadratic function to the background in the immediate vicinity of the 
peak, and subtracts this fitted background to obtain the peak sum. 

Branching ratios were determined by first fitting these peak 


intensities to the function 
W(8) = ag(1 + agP.(cos 6) + a,P,(cos 6)) (3-1) 


taking into account absorbtion of the y rays in the 0.005'' Ta backing 
and the 0.055'' Al plate, which were at an angle of 30° with respect 
to the beam axis. The ag coefficient of this expansion was then 
corrected for the efficiency of the Ge(Li) detector, which was 
measured using °’Co, *2Na, !37Cs, °*Mn and ®°Co sources of known 
intensity. The efficiency curve for the 48 cc detector is shown 
in Figure 9. 

The population parameters of the initial state in the y decay 
were calculated from the Statistical Compound Nuclear computer code 
AK (Gr 71b) which is based on the code MANDY of Sheldon and 


Strang (Sh 69), but expanded to include calculation of the y-ray 


ae A a 
to. yersns cay 6 Pre <tottetanio yore orld a 
Youd (86 9713 3 1w esansbionios fea efT  .722 to shen. as bne hie te 
nee ytteeto sd neo 2noisiensit Ved etd bis Vex 830 , vod 13d Ved < aa 
sua i a 
evBtqeib tasuqmos sda mort yisoozib benisido stew esisiaensini ysr-y : i 
30 Asani l ‘6 tadidia 2217 ote’, (ANY. 70) THMMU2 mereong: ofa gpiee 

ahs 70 ysiniaiv sisibenm! art nl bnuotpdoed oft oF notion? aigerbeyp 
muz dee sft nistdo of bruorpioed besti? 2idt eioeizdee bam a 
isbq s2od3 pnittit +2117 yd benimiesab syew eoftes pnifonevs : 
noljonut sdz 07 zeizienaant " 


e) ((8 209) 0918 + (6 2ande%gs * Hye = (OW i 


oy 

pniaaed sT "800.0 sft ni eye, Y ed? Yo aolsdioeds touoadss ont palates ~ 
#2sq261 dtiw “Of to sipre ne 36 oisw dairw .esefq fA "240.0 edt wer 

ners 2ew noleneqxo zing To tnaisitiaoa 96 sit .2ixe meed ot oF 2 


26w doinw .rolsetsb (11)s9 edt to yonsiaitts sd? 107 betoe1109 it 
eer 


nwond Yo esa7u0e 03°? bne aM*® .29°F4 .cWS* ,o0** eniew beruesem 
at 

nworle 21 s0399%5b 29 8) edz 307 gvius yonsisitie aft .yrtenesnt 
ee Si 


(yeaah bgt i ni stese Tetzini ena to etasame7eq he) aes - 
‘44 oo” 
sto2 1s2uqnoo testsu# bavogmo3 Jeoiteitere ada mon? bezstuolss 2: , 


bne nobiad2 to Vowan sboo. ‘orld no baeed 2i dotdw (APY 13) HA” 
yey ans Fo. eal stort er tud de aan 


600 


159>0 


RANDOM COINCIDENCES 


400 


200 
0 
200 400 600 800 
ENERGY (keV) 
5 
600 un 


REAL + RANDOM COINC! DENCES 


66 3+278 


400 


200 


200 400 600 800 
ENERGY (keV) 


Figure 8. A typical y - y coincidence spectrum, showing y rays 
detected in the Ge(Li) which are in coincidence with 
the 278 keV transition in ©*Cu. This spectrum was 
obtained for the Al geometry, at 6=55°. 


34 


(Vea) ‘YoASHa 


239W39019K103 MOGHAA + JAA 


eess\9 


102 ; 
| 10 
> 
oO 
Ps, 
Lid 
Oo 
i 
= 
Lid 
ce 
ke 
<x 
—_ 
Lu 
ce 
lo~ 
Ome | | gem ei 
0 On? 0.4 0.6 0.8 120 Vee ina 
y RAY ENERGY (Mev) ————~» 
Figure. 9. 


Efficiency curve for the 48 cc Ge(Li), measured with 
57C€o,22Na,137Cs,54Mn and ©9Co sources. 
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linear polarization coefficients as well as angular distribution 
coefficients. The transmission coefficients required as input to 
this program were calculated by the computer program ABACUS (Au 64), 
using the optical model parameters of Bjorkland and Fernbach (Bj 58) 
for both proton and neutron channels. Partial waves up to and 
including 2 = 4 were considered in all channels. 

For many of the transitions in °*Cu, both angular distribution 
and y - Y triple correlation measurements in two geometries were 
performed. For these transitions, all three sets of data were 
used to minimize a total x? to find acceptable values for the spin 
and mixing ratio, using the computer code SUPERCHI (Gr 73). Since 
only three linear polarization measurements were con ehey were not 
included in the computation of the total y* for the respective decays, 


and they were analyzed separately. 
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CHAPTER 4 
EXPERIMENTAL RESULTS 


A typical Ge(Li) singles spectrum, taken at a proton bombarding 
energy of 3.75 MeV, is shown in Figure 2. y-ray transitions between 
levels in ®*Cu are identified by the energies of the initial and 
final states in the decay. The energies determined in this work agree 
well with the more accurate results of Shera and Bolotin (Sh 68), and 
since the emphasis in this work was not on precise energy determination, 
the values of this reference have been adopted throughout. 

Figure 10 summarizes the decay scheme of the levels below | MeV 
excitation in ©*Cu, as well as the branching ratios and spin-parity 
assignments determined in this work. This decay scheme differs from 
that of Shera and Bolotin (Sh 68) in the following respects: 

i) A weak 362 keV y ray is observed, which is assigned 
to the ground state decay of the 362 keV level. 
ii) A 719 keV y ray is observed, which is attributed to 

the 878 + 159 transition. 

iii) The 534 keV transition, which Shera and Bolotin (Sh 68) 

assign to the 895 > 362 decay, is re-assigned as the 

878 + 344 transition, on the basis of coincidence with the 

344 keV y ray. 

iv) A 736 keV transition is observed, which is attributed 


to the 895 > 159 decay. 
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Figure 10. Level scheme for ©*Cu, showing spin-parity assign- 
ments and y-decay branching ratios determined in 
this experiment. 


4.1 Branching Ratios 


A summary of the branching ratios measured in this experiment 
is given in Table Il. Although several y-ray studies of ®*Cu 
have appeared in the literature, only Shera and Bolotin (Sh 68) give 
results for the branching ratios. Their results, as well as those of 
Drum (Dr 70) who studied only the levels up to 663 keV, are also shown 
for comparison. The results of Shera and Bolotin (Sh 68) have been 
corrected for their improper placement of the 534 keV decay. 

For some of the weaker decays in ©*Cu, angular distribution 
measurements could not be made, so it was not possible to use the 
procedure outlined in Section 3.5 to calculate the branching ratios. 
For these transitions, it was assumed that the yield at 6 = 55° gave 
a good estimate of the a, coefficient in equation (=|) This will 
be the case provided that the a, coefficient is small, since 
P, (g=55°) = 0.0. 

The three sets of results in Table I! are in good agreement, with 
the exception of the 895 and 927 keV levels. For the 895 keV level, 
Shera and Bolotin do not observe the 736 keV transition to the first 
excited state, which we measure to be of comparable strength to the 
branch to the 278 keV level. This decay was also observed by 
Davidson et al. (Da 70), but they give no estimate of the relative 
strength of this branch. 

For the 927 keV level, our branching ratios are only approximate 
values, and it is difficult to assign reliable error estimates to the 
values quoted. The uncertainty regarding this level stems mainly 


from the presence of the strong 770 keV y ray from the °*Ni(p,yy) ©°Cu 
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159 
278 
344 


362 


574 
609 


663 


139 


746 
878 


895 


927 


A summary of the y-decay branching ratios for levels below 
1} MeV in ©4Cu. 
compared with those of Shera and Bolotin (Sh 68) and 
Drum (Dr 70). 
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reactions which partially obscures the 768 keV transition from the 

927 to 159 keV level. Estimates of the contribution of the 768 keV 

y ray to this unresolved doublet were made by fitting two Gaussian 
peaks to the spectrum at 6 = 55° using the computer code PEAKFIT 

(Ea 72), but the overall fit to the peak shape was poor. A similar 
procedure was followed for the 927 keV to ground state transition. 

In Figure 2, the peak labelled 927 + 0 is about twice the width of the 
other ©*Cu peaks, and it is likely composed of two y rays. The 
source of the other transition has not been identified. Because of 
the above uncertainties, the branching ratios of Shera and Bolotin 


(Sh 68) for this level are considered more reliable. 


i2 The Levels of °*Cu 


Table I11 gives the experimental a, and a, coefficients for the 
measured angular distributions of transitions between low-lying states 
in ®*Cu. Table IV lists the measured a> and a, coefficients for 
both geometries of the y - y triple correlations. The triple 
correlation results were obtained at a bombarding energy of 3.75 MeV. 

A summary of the spins and mixing ratios determined in the 
present work is given in Table V. Table VI compares these spin 
assignments with those obtained in previous works, and Table VI] 
compares the mixing ratios measured in this work with those of 
Davidson et al (Da 70). 

lt was mentioned earlier that, in many cases, more than one 
measurement was made at each ancien Although these measurements 


were treated as independent data points in the y* analysis, the 


TIFHARS sboo eae sit enizu "de.= 0 te. accel a os é 4, 
aa yh. 
wstimi2 A .100q 2sw sqere Assq Say oF 41? flevevo stz gud hata 
.moigiens1t aiész2 bnuotp oF Ved {SC sf? 107 bewollot 2ew o7ube207q 
aia 2o dtbiw od¢ soiws tuode 2] 0 + {SC befistsl wseq ay <S ewweit at aA 
sat .e¥e1 yowt To baeogtios yledif 2t t+? bre ,2teeq u3** vetso = 


$9 s2u6089 .beitigasbi need ton 260 soitlensss ved¥o ada to sowoe 


nitoltod bie 619c¢ 4o 2oite: onidariend sno ,asisnietisonv svete it 


-stdeiier siom bersbienco eis Isvel aida 30? (83 az) 7 
f 
: 


ui*? to al evos aT S.A 


sti] 10? 2tasizltisos 4s bose .6 (stasmizeqxs sit 2evig Ii! sideT 


asset pity! -wo!l naswtod enoitienst3 to enolsudivzeib seluens bsiwesom - 


ot etrsiciiteoa ,& bas 6 beweeem eft esei! VI stdetT .wd*? at 
glqiid sit  .znolseierios siqitt y - y 93 to esiviemose riod 
.VaM 2¥.£ Yo yorene pdibiedmod 6 16 banistdo stew etfuesy not satene _ 
sia al banimtvezsb zoltat poaixim base eaiqe ada to yremme A 
niq2 seont zs1eqmoo [V sfdeT -V sidsT ni nevip at stow inses1q ay 
(1V sfdeT bos ,2xnow euoiveig ak banistde eens daiw zinamngizes 7 
to szort dsiw stow cirls ni bawwessm 2oltey onixim oft 2s18qmo> 
(0% 00) {e629 noebivad ” 
sno ned 210m . 28265 ynsm ni , ied? yellieo bsnoisnem 2ew 3! 7 
ztnomevuessi seedy dguodsziA stone i369 16 sbsm 26w tnemstu: | 
ars -etayleps Ay. ait ni etniog 466 Jnebheqabni 28 bs | 


Table I!! Experimental az and a, coefficients for the y ray 
angular distributions. 
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weighted average of all measurements at one angle is plotted in the 
diagrams. 

Both the angular distribution and angular correlation measurements 
are independent of the parity of the rade kion. Since they are not 
sensitive to the direction of the electric vector of the radiation. 
These experiments may be sensitive to the parity of the initial state 
in the decay, however, if the degree of alignment (specified by the 
BL (da) coefficients) is parity dependent. This is the case if the 
Statistical Compound Nuclear model is used, since the population 
parameters are determined, in part, by the optical model transmission 
coefficients 193 which depend on the orbital angular momentum of the 
incoming and outgoing particles. 

Figure 11 shows this dependence on parity for the angular 
distribution of the 278 keV transition. The effect is most marked 
for low spins (J = 1 in this case, since J = 0 always gives an 
isotropic angular distribution). For all of the transitions studied 
in this work, calculations were done for both parity assignments, and 
in no case was it possible to exclude one parity assignment on this 
basis. Since the differences are usually small, al] results shown in 
this section (except for the linear polarization measurements) are for 
positive parity assignments. The measurements of Park and 
Daehnick (Pa 69) indicate that there are no negative parity levels 


in ©*Cu below 1.4 MeV excitation. 
4.2.1 The 159 and 278 keV Levels 


Both of these levels decay 100% to the ground state, which has 


spin-parity J" = iy (Le 54, Do 66, Wa 67). The measured angular 
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Figure Il. 


tan! § 


x? curves for the 2780 angular distribution, 


showing the dependence of the angular distribution 
on the parity of the initial state in the decay. 
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Figure 12. 


tan! § 


The 159 keV level: a) x2 curves and b) angular 
distribution for the 1590 transition. 
The 278 keV level: c) x2 curves and d) angular 
distribution for the 278+0 transition. 
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Figure 13. 


The 278 keV level: a) .x% curves and b). linear 
polarization distribution for the 278. keV y ray. 
The solid lines in b) show the expected linear 
polarization distributions for values of J. and 
6 which gave good fits to the angular distribution. 
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distributions and y* curves for these two transitions, shown in 
Figure 12, indicate a unique spin assignment of J = 2 in both cases, 
with a mixing ratio of §6 = 0.12 + 0.04 for the 159 keV transition, 
and § = 0.10 + 0.02 for the 278 keV transition. 

Figure 13 shows the results of the linear polarization 
measurement for the 278 keV transition. The solid lines in the 
lower half of this diagram are the expected linear polarization 
distributions for the various spin and mixing ratio combinations 
which give a good fit to the angular distribution data. The 
J" = oe 6 = 2.14 assignment is considered unlikely because of the 


large M2/E1 enhancement required to fit the data, and we therefore 


assign positive parity for this level. 


4.2.2 The 344 keV Level 

This level decays primarily to the ground state, with a 4% 
branch to the 159 keV level. The angular distribution of the ground 
state decay is isotropic within experimental limits (a2 = 0.01+0.03) 
and acceptable fits were obtained for J = 0, 1 and 2 (Figure 15). 

The measured angular distribution of the 185 keV transition was 
also isotropic, but the Al geometry of the 344 + 159 + 0 angular 
correlation had an appreciable anisotropy of ap = -0.17+0.06. The 
curves shown in Figure 14 indicate acceptable fits for J = 1, 2 or 3, 
but J = 3 is not allowed from the ground state angular distribution. 

The J = 0 assignment can also be eliminated from.the 609 + 344 
angular distribution (Figure 19). The spin of the 609 ey) bevel is 


known to be 2 from analysis of the ground state angular distribution 
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The 344 keV leyel: a) x? curyes, b) angular distri- 
bution for the 185 keV y ray, c) Al. geometry. and 
d) C1 geometry of the 344+159%0 angular correlation. 
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The 344 keV level: .a) x2 curves and b) angular 
distribution for the 344+0 transition. 

The 362. keV level: c) x2 curves and d) angular 
distribution for the 362+159 transition. 
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Figure 16. 
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The 344 keV level: a) x* curves and b) linear 
polarization distribution for the 344 keV y ray. 
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(Section 4.2.5), and the angular distribution of the 265 keV transition 
is inconsistent with a J = 0 assignment for ene 344 keV level. 

The J = 2 assignment can be ruled out on the basis of a 
combination of the angular distribution and linear polarization 
distribution of the 344 keV y ray. For J = 2, an isotropic angular 
distribution occurs for a mixing ratio of 6 = -0.19, which gives rise 
to a polarization of about 30% at 6 = 90°. The near isotropy of both 
the angular distribution and the linear polarization distribution 
(Figure 16) is consistent only with a spin-parity combination of 


jl aes 


4.2.3 The 362 keV Level 

This level decays 98% to the 159 keV state, with a weak ground 
state branch also being observed. The angular distribution of the 
203 keV y ray, shown in the lower half of Figure 15, results in a 


unique spin assignment of J = 3, with a mixing ratio of 6 = 0.06+0.03. 


4.2.4 The 574 keV Level 

This level is only weakly excited in the (p,ny) reaction, and 
was observed to decay only to the J = 3 state at 362 keV. Although 
the presence of this state has been reported by many authors, no 
previous angular distribution measurements have been reported for 
the 212 keV y ray. 

In the present study, both the angular distribution and the 
574 + 362 + 159 angular correlation were measured. The results of 
these measurements are presented in Figure 17, showing possible fits 


for J = 4 (§ = 0.0140.03) and J = 2 (-4.7 < Of fe. G)) 
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Figure 17. 


The 574. keV level: .a) x? curves, b) angular 
distribution for the 212 keV y ray, c) Al 
geometry and.d) C1 geometry of the 574+362+159 
angular correlation. 
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Table V A summary of the spin assignments and multipole mixing ratios 
determined in this work. 


(1) 
EF E J 6 
159 159 2 Ofl2e= Or os 
278 278 2 (+) 0.10 + 0.02 
344 344 1s All values 
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362 203 3 0706812 0.02 
57h 212 h O70 let 10203 
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609 609 2 0.30 + 0.08 
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Table VI A comparison between the spin assignments determined in the 
present experiment and those of previous works. 
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The J = 4 assignment seems more likely, for the following reasons. 


The level is very weakly populated in this reaction, suggest ing that 
the spin is high, and that particles of high angular momentum (and 
hence smaller transmission coefficients) are required to form this 
state. The fact that the only observed y decay is to a J = 3 state 
also tends to favour the assignment of J = 4. The §&Zn(d,q) 
measurements of Park and Daehnick (Pa 69) indicate L = 4 for this 


level, limiting the spin to a 4 or i. 


4.2.5 The 609 keV Level 

The primary mode of decay for the 609 keV level is the ground 
State decay (82%) with weaker branches of 8%, 4%, and 6% respectively 
to the first three excited states. The angular distribution of the 
ground state branch was highly anisotropic (a, = -0.34+0.01), and 
the y* analysis resulted in a unique assignment of J = 2 (Figure 19), 
with possible values for the mixing ratio of § = 0.30+0.08 and 
§ = 1-1520.25. 

The linear polarization distribution for this decay is shown 
in Figure 20. Of the four possible solutions for the fit to the 
angular distribution (two possible mixing ratios for each parity), 
the 274 § = 1.15 and oan § = 0.25 solutions can be eliminated on 
the basis of the polarization measurement. The J" = ee 6 = 1.30 
assignment is considered unlikely because of the large M2/E] 
enhancement required to fit the data, and positive parity is the 
preferred assignment. 

The measured angular distribuiton and y - y angular correlation 


of the 450 keV transition are shown in Figure 18. Here again, the 
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Figure 18. 
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The 609 keV level: .a) x* curves, b) angular 
distribution for the 450 keV transition, c) Al 
geometry and.d) C1 geometry of the 609°1590 
angular correlation. 
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Figure 19. 


The 609 keV level: a) x2 curves and b) angular 
distribution for the 609 keV y ray.. 

c) x2 curves and d) angular distribution for the 
609+344 transition. For c) and d), J refers to 
the spin of the 344 keV level. 
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The 609 keV level: a) x* curves and b) linear 
polarization distribution for the 609 keV y ray. 
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J = 2 assignment fits the data best, with the J = 3 assignment being 
almost eliminated on the basis of the 0.1% confidence limit. The 
mixing ratio of the 450 keV y ray is found to the 0.10220207. 

Angular distribution measurements were not possible for the 


331 keV transition. 


4.2.6 The 663 keV Level 

This level was observed to decay almost equally to the ground 
state and first two excited states, with a weaker branch (6%) to the 
344 keV level. | 

The angular distribution of the ground state branch was isotropic 
(a, = 0.00+0.02), and as was the case for the 344 keV transition, 
acceptable fits were found for J = 0, 1 and 2 (see Figure 22). For 
the J = 2 assignment, the values of the mixing ratio which resulted in 
isotropy were quite limited. Figure 22 also shows the measured 
angular distribution of the 319 keV y ray (a. = -0.02+0.02), which 
likewise gave possible fits for J = 0, 1 and 2. Measurements on the 
504 keV transition were not possible, due to the proximity of the large 
annihilation radiation peak at 511 keV. 

The angular distribution and two geometries of the angular 
correlation for the 385 keV y ray are shown in Figure 21. Both the 
distribution and the Cl geometry are isotropic, indicating that the 
initial state is not strongly aligned, while the Al geometry has 
ap = -0.18+0.03. Simultaneous analysis of these measurements yields 
a unique assignment of J = 1 for this fevers with § = 0.07+0.05 or 
oH] 5 6 5 72.99) Re Ihis assignment is in disagreement with Park and 


Daehnick (Pa 69), who favour A 30 assignment on the basis of an L = 4 
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Figure 21. 


The 663 keV level: a) x curves, b) angular 
distribution for the 385 keV transition, c) Al 
geometry and d) C1 geometry for the 663+278+0 
angular correlation. 
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Figure 22. 
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The 663 keV level: a) x2 curyes and b) angular 
distribution for the 663 keY transition. 

c). x2 curves and d) angular distribution for 
the 663+344 transition. 


Table VII A comparison between the mixing ratios determined in the 
present experiment and those of Davidson et al (Da 70). 


E § 
Spins Present Work Da 70 
159 27>] 0.12 + 0.04 -0.01 + 0.02 
278 2> 1) 0.10 + 0.02 0.025 + 0.025 
3hh 17] All Values (C0) 
185 172 OnrOtUr Ouse lcocas st - 
203 3>2 0.06 + 0.03 G04 02035 
609 27> 1) 0.30 + 0.08 On0 see Un 0s 
450 2>2 0.02 + 0.07 -0.05 + 0.10 
265 2> 1) 0.24 + 0.17 
663 1> ] All values = 
385 1>2 G07 0205p 0 -- 2.9 ~ 
Salle hues One xe onsen Gy . 
580 2>2 -0.18 + 0.11 0.05 + 0.04 
3>2 -0.38 + 0.04 -0.34 + 0.06 
46] 2>2 -0.29 + 0.25 Oe ua a1 
342 -0.43 + 0.10 =0.15 = 0.07 
377 2iPs 3 -0.11 + 0.18 ~ 0.05 
Ses 0.57 + 0.18 =0..50 
468 3>2 0.08 + 0.03 -0.04 + 0.035 
736 0, b0e2, 0.13 
617 3>2 0 <O7cegEs ty 2015 OL Ow +2 Ole 
649 1>2 ONOSFOS TIN =5: /56<-2.5 ~ 


S02 i005 


aso.0 + 20.0 
(o < 8) 


2£0.0 = 40.0 


£010 +18 
eee ee 


=r rane = ei > 
: Picks a La, 


0.0 + S10 fs, : eet 
$0.0 + 05.0 f+s a) 
feV [TA oe | aee 
O.g-2821.0- Of 0208.0 Nan . 
£0,0 + 30,0 s+ef £08. 
wee 1+s ~ © 0a 
0.0 + 86.0 S+§ oe 
ttt ds.0 “2 | co a 
ins 
s y a 0. . 7a 
> a , bet | ia o 
a 7 
ine ae oe i 08 i 
0.0 + al~ * ee a 
Bios eh.ee on ST Be 
f.0-= % + - 7 
| aos a it~ Paul we 
20+ V2.0 +t ) Ga 
= 8 fl 
£0.0 + 80.0 S+e£ if el 
.0 + 08 Bey 
> oe, g+f€ 
2. 84e3st.2- iad g*T 
= ay 
uf Pa _ a 
y — . » 


64 


transfer from their ©®Zn(d,a) measurement. Wellborn et 'al (We 71) 
studied the 663 + 278 + 0 cascade, and found acceptable fits for either 
J = 1or J = 2 for this state, with the J = 1 assignment more probable. 


Bass and Stelson (Ba 70) also foundJ = 1 (2) for this level. 


4.2.7 The 739 and 746 keV Levels 

The 739 keV level decays 77% to the first excited state, 6% 
to the 278 keV level, 4% to the 344 keV level and 13% to the 362 keV 
level. No ground state branch was seen. 

The results for the 580 keV decay are shown in Figure 23, and 
indicate either J =27) dea eaOielGsOe11, or k= 3, 6 = -0.3830,04., 
Similarly, the analysis of the results for the 461 keV y ray yield 
possible solutions for either J = 2 or J = 3 (Figure 24), as does the 
angular distribuiton of the 377 keV transition, shown in Figure 25. 

For the 746 keV level, only a single branch to the 278 keV 
state was observed. The results of the measurements on the 468 keV 
transition are shown in Figure 26, and the x* analysis results in a 


unique assignment of J = 3, with a mixing ratio of 6 = 0.08+0.03. 


4.2.8 The 878 keV Level 

This state was observed to decay 55% to the ground state and 
42% to the 344 keV state, with a very weak (3%) branch to the 159 keV 
level. The absence of strong decays to the two low-lying J = 2 levels 
suggests J = 0 for this level. 

The angular distributions of the 878 keV and 719 keV y rays are 
shown in Figure 27. Both of these decays were almost isotropic 


within experimental error (a, = -0.05+0.04 and a, = -0.16+0.15 
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Figure 23. The 739 keV level: a) y* curves, b) angular 


distribution for the 580 keV y ray, c) Al 
geometry and d) Cl geometry of the 739+1590 
angular correlation. 
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The 739 keV level: a) x* curves,’b) angular 
distribution for the 461 keV y ray, c) Al 
geometry and. d) C1 geometry of the 739+278+0 
angular correlation. 
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Pigure 25. The 739 keV level: a) x2 curves and b)--angular 
distribution for the 739+362 transition. 
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Figure 26. The 746 keV level: a) x2 curyes, b) angular 
distribution for the 468 key y ray, c) Al 
geometry and.d) C1 geometry of the 746+278+0 


angular correlation. 


respectively), and the x* analysis shows acceptable fits for both 
distributions for J = 0, ole ose 

The results for the 534 keV y ray are shown in Figure 28. The 
Al geometry of the angular correlation showed appreciable anisotropy 
(a, = 0.27+0.06), and the data result in a unique assignment of J = 0 
for this state. 

The J = 0 assignment is in agreement with the most likely 
assignments of Park and Daehnick (Pa 69) and Bass and Stelson (Ba 70), 


but not with Davidson et'al (Da 70), who find J = 1 or 2, with J = ] 


favoured. 


4.2.9 The 895 and 927 keV Levels 

The level at 895 keV was found to decay almost equally to the 
159 keV and 278 keV states (40% and 49% respectively), with an 112% 
branch to the ground state as well. Angular distribution measurements 
were performed for the two strongest branches, and are shown in 
Figure 29. Both measurements resulted in a unique spin assignment 
of 3 for the 895 keV state, with mixing ratios of 6 = 0.404+0.13 for 
the 736 keV transition, and 0.07 < 6 < 2.5 for the 617 keV y ray. 

As was mentioned in Section 4.1, accurate branching ratios for 
the 927 keV level could not be determined from our experiment. The 
strongest branch was to the 278 keV level, with weaker decays to the 
ground state and first excited state also being observed. 

The results of the angular distribution and angular correlation 
measurements for the 649 keV transition are shown in ei gure 30. 
Possible fits were obtained for J = Oee. 020420. 11 gon 5.7.5) 6 < ~2.5, 


| ; 
and J = 3, § = -0.114+0.05. Park and Daehnick (Pa 69) assigned J” = ] 
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Figure 27. 


The 878 keV level: a) x* curves, b) angular 
distribution for the 8780 transition. 
c)- x?! curves and d) angular distribution for 
the 878159 transition. 
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Figure 28, 


The 878 keV level: a) x2 curyes, b) angular 
distribution for the 534 keV y ray, c) Al 
geometry and d) CJ geometry of the 878+344+0 
angular correlation. 
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Figure 29. The 895 keV level: a) x2 curves and b) angular 
distribution for the 895>%278 transition. 
c) x* curves and d) angular distributlon for the 
895+159 transition. 
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Figure 30. The 927 keV level: a).x* curves, b) angular 
distribution for the 649 keV y ray, c) Al 
geometry-and.d) C1 geometry of the 927+278-0 
angular correlation. 
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for this level, on.the basis of a strong L = 0 transfer in the 
®67n(d,a)°*Cu reaction while Bass and Stelson (Ba 70) find J = 2 
most likely, with 0 or | assignments also possible. The J = 1 
assignment therefore seems more probable, although our data cannot 


convincingly distinguish between the two possibilities. 
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CHAPTER 5 
THE INTERMEDIATE COUPLING MODEL 


5.1 Motivation 


There is evidence, from the (d,p) stripping measurements of 
Park and Daehnick (Pa 69), for considerable fragmentation of the 
neutron single particle orbitals in ®*Cu. Although most of the 
Faro strength has been assigned to the two levels at 362 keV and 
574 keV excitation, considerable ie = 3 admixtures are also observed 
for the ground state, first two excited states, and the level at 
895 keV excitation. Much smaller fragments are observed for several 
States between 1.5 MeV and 3 MeV. 

The le = ] strength is concentrated mainly in the states below 
1.5 MeV, but is split among about ten energy levels in this region. 


As in the case of the f orbital, smaller fragments are observed 


5/2 
up to 3 MeV. 

As mentioned in Chapter 1, the low-lying energy levels of the 
even-even nucleus °*Ni exhibit a characteristic vibrational spectrum, 
with a 2° first excited state lying midway between the oO” ground 


he a2eande4 |aeSuch nuclei 


State and a triplet of states with J" = 0 
can be described by the collective vibrational model. The 
measured B(E2) for the 2° + 0° transition in ©2Ni of 154 + 17 


e*fm* (Du 66) is about 1] times larger than a single particle estimate 
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would give, again indicating that some sort of collective motion 

is taking place. anaes Gomez (Go 71) and Markham (Ma 71) have 
presented calculations in which the odd-A nucleus S4Cu is described 

as a single proton coupled to the surface vibrations of ®2Ni. Both 
of these calculations are able to reproduce the energy spectrum of 

®3Cu extremely well, as well as the single particle spectroscopic 
factors and B(E2) values for the y decays. In view of the success 

of these calculations, it is tempting to apply a similar model to ©*Cu. 

Calculations based on the Intermediate Coupling model have been 
carried out by several authors, mostly for odd-A nuclei (Ra 63, Ch 67, 
Ca 71, Go 71, Ma 71). Of the calculations which have been performed 
for even-A nuclei, the majority have been concerned with describing 
the properties of even-even nuclei in terms of two identical particles 
coupled to the vibrating core (Ma 63, He 67, Pa 71), and only the 
calculation of Heyde and Brussaard (He 68) for ?*°La attempts to 
apply the model to odd-odd nuclei. This calculation gave fair 
agreement with the experimental energy spectrum of !*°La. 

Section 2 of this chapter gives a brief discussion of two 
possible approaches to calculations of this type, the ''Thankappen and 
True'' model, and the Intermediate Coupling model, and a discussion of 
the relative merits of each. Section 3 contains a detailed discussion 
af the Intermediate Coupling model for odd-A and odd-odd nuclei, and 
Section 4 gives the results of the calculations for ®*Cu. A 


discussion of the results is presented in Section 5. 
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5.2 Basis for the Model 


There are two general methods of approaching calculations in 
the Unified model. The first considers a single nucleon coupled 
to some sort of core erates but makes no assumpt ions about either 
the nature of these core states or the details of the interaction 
between the core and the extra particle. This type of approach 
was suggested by Lawson and Uretsky (La 57) and de Shalit (Sh 61), 
who assumed only that the particle-core interaction energy could be 
expressed as a sum of scalar products of two tensors, one of which 
acts only on the co-ordinates of the particle and one of which acts 
only on the co-ordinates of the core 

thie 3 T (p) *T, (e) (5-1) 
This model was later developed by Thankappen and True (Th 65, Th 66), 


who used an interaction of the form 


Hive = ~E((c)*F(p)) - n(Q(c)-Q(p)) (5-2) 


> 
where J is the total angular momentum operator and Q the quadrupole 


moment operator. Using this model, they calculated the properties 
of the low-lying states of ®3Cu and obtained good agreement with the 
experimental energy spectrum and B(E2) values for this nucleus. 

In the Thankappen and True model, the matrix elements of the 
total Hamiltonian depend on the interaction strengths n and &€, and 
reduced matrix Bvenente of the form <j'||Q(p)|[j> and <J'|/Q(c) | [J>. 


The particle reduced matrix elements <j'||Q(p)||j> are calculated 
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using harmonic oscillator wave functions. Since the model does not 
assume a detailed knowledge of the core erates it is not possible to 
calculate the core reduced matrix mienenten and they must be treated 
as adjustable parameters. One obvious disadvantage of this type of 
approach is that it is usually necessary to restrict the number of 
core states in the calculation rather severely, in order to keep the 
number of free parameters reasonably small. 

The second approach to calculations of this type is to assume 
that the core is capable of performing collective oscillations, of 
the type described by Bohr and Mottelson (Bo 60). Usually, it is 
assumed that the core is described by a quadrupole surface 
oscillation, and that the form of the interaction between the extra 
particle and the core is known. The only adjustable parameter in 
this type of calculation is the strength of this interaction. 
Calculations of this type have been used by several authors (Go 71, 
Ma 71) to describe odd-mass nuclei in the vicinity of the Cu isotopes. 

We will follow the second approach in studying the odd-odd 
nucleus ©4#Cu. Besides the obvious computational advantage of having 
fewer adjustable parameters in the mode], it is felt that the added 
generality of treating the ®*Ni core states as unknown is unnecessary. 
The spectrum of ©*Ni, along with the success of the calculations 
fOhus OCU, seems sto ind! cate) that the °2Ni core is adequately described 


in terms of the vibrational model. 
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5.3 Formalationsef the Model 

This section contains a brief description of the mathematical 
formulation of the Intermediate Coupling model. A more detailed 
discussion of the model, including the derivation of most of the 
formulae presented here, is given in Appendix B. 

For simplicity, we first consider the case of a single particle 
coupled to a core. It is assumed that the core is capable of 
performing collective oscillations, and that the nuclear radius can 
be expanded in the form 

R(6,6) = Roll + : ay Y," (8,6) ] 1S) 

u 

The term with A = 0 corresponds to changes in the density of the 
nucleus. Such ''breathing modes'' of excitation require a great deal 
of energy, and are not expected to be important in the description 
of low-energy properties of the nucleus. The A = 1 term corresponds 
to oscillations of the centre-of-mass of the nucleus. The first 
term which is expected to be important for the low energy properties 
of nuclei is that for A = 2, and we will assume that all higher order 
terms in (5-3) are negligible. The extra particle is assumed to have 
several single particle orbits available to it. 

The total Hamiltonian for this system can be written as a sum 


of three terms 


H=H + H, + Ais (5-4) 
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The Hamiltonian for the core (Ho) can be Written as 
H = (bh -b -+1/2) tw (5-5) 
Cae HU 


where OT and Ze are respectively the creation and annihilation 
operators for (spin 2) phonons of energy tw. The explicit form of 
the Hamiltonian for fhe extra particle (H.) is not important for 
this discussion, and it is sufficient to note that when it acts a 

a single particle wavefunction, it has the single particle energy a 
for that particular orbital as its eigenvalue. The interaction 


Hamiltonian is given by 


ls at rv 7 
He = (=) hwé [ (-) b+ b1Y2 (9,6) (5-6) 


where (8,6) are the polar angles of the extra particle. The 
strength of this interaction is given by the dimensionless coupling 
constant &. 

The basis states are chosen as those of the uncoupled system, 
and are specified by the number of phonons in the core state N, which 
are coupled to total angular momentum R, the single particle angular 

> 
momentum j and the total angular momentum | = [R a Tle Neglecting 
the zero-point energy, the diagonal elements of the total 


Hamiltonian are 
<NRJ31M[H. + H INR 5 1M = Nw + e, (5-7) 


The off-diagonal elements in the Hamiltonian matrix are due entirely 


to Hate and including this term results in mixing between the basis 
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states. The selection rules for non-vanishing matrix elements of 
Hint arerun = | eARecre wn les 2 and) AL=10 5 2 only. The eigenvectors 
of the complete system are given as linear combinations of the basis 


States 


JE;IM> = 5 BEY 
NRj 


NRJ1) | NRJ 3 1M> (5-8) 


and the expansion coefficients B™ (NRj 1) are found by diagonalizing the 
Hamiltonian matrix. 
For a system composed of two particles outside a collective core, 


the Hamiltonian is 


Hatley H (1) + Hagel, + H. (2) see Nig (5-9) 


where the indices 1 and 2 refer to particles ] and 2 respectively, and 
H,, is the residual interaction between the two extra particles. As 

explained in Appendix B, it is convenient to choose as a basis set 

for this system, states in which particle 2 is coupled to eigenstates 


of the core plus the first particle 
|a,Jij2;IM> = X (JajasM-me mo|1M) |oyJaM-m2 >| j2m2> (5-10) 


The states |v. J yM-m2> are given by equation (5-8), and are found by 
diagonalizing the first three terms in the Hamiltonian (5-9). The 


first four terms in (5-9) are diagonal in this representation 
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where ej is the single particle energy of particle 2, and es is 
2 1 
the energy of the core plus first particle eigenstate. H sayy (2) and 


Hi2 have both diagonal and off-diagonal matrix elements in this 
representation. Explicit formulae are given in Appendix 8. 
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5.4 Application of the Model to °‘*Cu 
5.4.1 Energy Levels of °?Cu 
The first step in this calculation is to consider the nucleus 


®3Cu as a single proton coupled to the ®2Ni core. This problem has 


been treated by Gomez (Go 71) and Markham (Ma 71). The odd proton 


is assumed to have the 2P 379° Ife72 and 2P 1/2 orbitals available. The 

i i i € & = |. V 
single particle energy spacings were taken as 5/2 3/2 1.500 Me 
and E12 - E372 = 2.200 MeV in agreement with the best-fit values 


found by Gomez. These values are not very different from those 
found by Blair (Bl 65) from the °*Ni(*He,d) reaction. |The phonon 
energy fw = 1.172 MeV was obtained from the excitation energy of 

the first ae level in the °7Ni core. In order to keep the number of 
components in the °°Cu eigenstates small enough so that the two 
particle catculations remain reasonably simple, only core states with 
N < 2 were Sere dare: in this calculation. 

Using the above values for the parameters, the Hamiltonian 
matrices for total spin values | = 1/2, 3/2, 5/2, 7/2 and 9/2 were 
diagonalized using the computer code CUPPLE-| (Gr 73a). Figure 31 
shows how the energy levels for the different spin values vary as a 
function of the coupling parameter &. The best fit was found for a 


value of € = 2.12 which is indicated by a vertical dashed line in 
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Figure 3], The Intermediate Coupling Model: Energy levels 


of the ©2Nit+tp system as a function of the 
coupling strength €. The best fit to the 
levels of &3Cu occurs for E=2.12. 
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Figure 32. A comparison between the experimental energy 
spectrum of ©&3Cu and the level scheme calculated 
from the Intermediate Coupling model. 
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Figure 31. Figure 32 shows the comparison between the calculated and 
experimental energy levels for ©%Cu, and Table VIII lists the expansion 


coefficients Be (NRj1) of equation (5-8) for the lowest energy levels. 


5.4.2 Energy Levels of ®*Cu 

For the ®4*Cu calculations, the lowest 11 eigenstates of the °%Cu 
system were used as basis states, and the odd neutron was assumed to 
have the lf 72 and 2P 179 orbitals available. The energy separation 
for these two single particle orbitals was chosen originally as 
E772 57 E5/2 = 0.500 MeV in agreement with the value deduced from the 
(d,p) measurements of Fulmer and McCarthy (Fu 63). This value was 
later reduced to 0.400 MeV to obtain better agreement with the 
experimental energy spectrum. The value of Eo (for the proton) was 
fixed at 2.12, the value which gave the best fit to the ®°Cu spectrum. 
The two-particle residual interaction H,, was chosen as a 6-function. 
The calculation then has three parameters, these being the neutron 
interaction strength ony the single particle energy spacing and the 
two particle interaction strength 

Figure 33 shows the calculated energy levels of ®*Cu, as a 
function of the neutron coupling strength ows FOG =) 0.0. The 
best fit corresponds to a value of on a oie No significant 
improvement in the fit was obtained when the two particle interaction 
strength n was varied from zero. The left-hand side of Figure 34 
shows the calculated energy levels of ©*Cu using this model, compared 
to the experimental spectrum shown in the centre. Only positive 
parity levels are shown, and spin assignments for the levels above 


1 MeV are taken from the work of Park and Daehnick (Pa 69). Table IX 
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The Intermediate Coupling model: Energy levels 
of the ©3Ccutn system as a function of the 
coupling strength €. The best fit to the 
levels of ©4Cu occurs for &=1.5. 
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Figure 34. A Co OR between the experimental level 
scheme of ©*Cu and the energy levels predicted 
by the [Intermediate Coupling model. 
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gives the expansion coefficients for the first few levels of °"Cu. 
Table X lists the best fit values for the parameters used in this 
calculation. 

An inspection of Figure 34 shows that not all of the levels of 
®"Cu are reproduced by this calculation. in particular, there are 
no levels in the theoretical spectrum which correspond to the 2" 
state at 278 keV, the J = 1 level at 663 keV or the J = 3 level at 
746 keV excitation. It was thought that these levels might be due 


to the excitation of a P3792 neutron into the f orbital, and could 


bYZ 


therefore be explained as a P372 neutron hole coupled to the levels 
of eesCu. 

The outline of this calculation follows that of the previous 
one. The parameters for °°Cu are identical to those used by 
Gomez (Go 71). Figure 35 shows the energy levels of e Cipas a 
function of ene with the best fit value of oS = 2.10 shown as a 
dashed line. Figure 36 shows a comparison of the experimental 
and calculated energy levels of “Cur 


The predicted energy levels of ®*cu arising from this 


configuration are shown on the right hand side of Figure 34. lt was 


not possible to obtain a good fit for any value of the coupling 
strength Sh) the levels shown in this diagram are those for 
cn Ed Alaein The agreement between the experimental levels and the 


theoretical predictions is not good. The parameters used for this 


calculation are given in Table X. 
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The Intermediate Coupling model: Energy levels 
of the S4Ni+p system as a function of the 
coupling strength €. The best fit to the 
levels of ©9Cu occurs for £2.10. 
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Figure 36. A comparison between the experimental energy 
spectrum of ©°Cu and the leyel scheme calculated 
from the [Intermediate Coupling model. 
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Table X Parameters used in the Intermediate 
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Coupling model calculations. 
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5.4.3 Spectroscopic Factors. 
Spectroscopic factors for the ®3Cu(d,p)®*Cu reaction are given 
by the expression 
Sa UR Meru ase -12 
bajy 7 UZ IAN3/2 505!) (5-12) 
where A°(3/2 j,;1) is the expansion coefficient for the ®*Cu 
eigenstate in terms of a neutron in the (2,,j,) orbit coupled to the 
J” = 3/2 ground state of ®3Cu. us is the ''emptiness'' of the orbit 
2 
VCore) forse Clinwandsis)eduall towl. for the Pisa orbital and 0.67 for 


the f orbital. 


Sy 


A comparison between the measured spectroscopic factors of 


Park and Daehnick (Pa 69) and the results predicted by the Intermediate 


Coupling model is shown in Table XI. The calculated values are in 
fair agreement with the experimental results, with the exception of 
the 278 keV and 609 keV levels. The values of Park and Daehnick 


have an estimated error of about 25%. 
5.5 Discussion of the Results 


Figure 34 shows a comparison between the experimental level 
scheme of ©*Cu and the level scheme predicted by the Intermediate 
Coupling model. The overall agreement for the levels below 1 MeV 
excitation is fairly good. The ground state and first excited 
State are well reproduced, although the excitation energy of the 
a level is about a factor of two too low, and the cy level at 
362 keV excitation is predicted very well by the model. The ia 


level at 344 keV is predicted to lie above the a but the energy 
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Spectroscopic factors (2J + 1)S for the ®3Cu(d,p) ©4Cu 
reaction, calculated from the Intermediate Coupling model. 
The calculated values are compared with the experimental 
results of Park and Daehnick (Pa 69). 


(D5 41) S 
J 
(Present Q Exp 
vera n (pa 69) THEORY 
] 1 0.17 0.11 
3 0.73 1.65 
2 1 0.84 Peete 
3 1.09 1.82 
2 1 0.85 - 
3 0.92 - 
1 1 1.18 1.92 
3 = 0.04 
3 ] 0.14 = 
3 Beg ey. 
h l = = 
3 3.64 5.18 
2 1 0.74 1.86 
3 - 1.00 
1 0.16 - 
4 ; Te 
(258) ] P's} 0.62 
3 3 - 0.00 
0 1 Only = 
3 a a 
3 ] 0.20 - 
3 0.89 0 
1 0.40 0.31 
3 - 0. 


These two levels were not resolved in the work of Park and 


Daehnick (Pa 69) 


+ Assuming J = 2 for the 739 keV level, and that this level 
corresponds to the J = 2 state predicted at 1134 keV. 
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discrepancy is only about 70 keV. The 2° at 609 keV and the oO” 
and zi at about 880 keV are also well reproduced by this calculation. 

This model fails to predict the second excited state, the 
J = 1 state at 663 keV, the doublet around 740 keV and the 927 keV 
level. Attempts to reproduce these levels by considering a P372 
neutron hole coupled to ®°Cu did not meet with success. These 
calculations do predict a closely spaced doublet with pio= ish rf 
which might correspond to the 663 and 746 keV levels, but the 2" 
member of this configuration lies about 250 keV above this doublet 
rather than 400 keV below it, which would be the case if this state 
corresponded to the 278 keV level. These calculations also predicted 
another low-lying O” State for which there is no candidate. 

The calculations could probably be improved by considering mixing 
between the particle and hole configurations. The level at 663 keV, 
which we have identified as being primarily due to a neutron hole 
configuration, decays with almost equal branches to the 159 and 278 keV 
levels, indicating that the configurations of these two low-lying 2* 
states are not radically different. The (d,p) stripping measurements 
of Park and Daehnick (Pa 69) give almost equal spectroscopic factors 
for these two levels, supporting the conclusion that these two states 
have a similar structure. The 895 keV level also decays to these two 
levels with comparable strengths. This is clearly not the complete 
picture, however, since it is then difficult to explain why the 
746 keV level decays only to the 278 keV level and not to the lower 
159 keV state or the ground state. This situation seems to suggest 


that this level is a relatively pure hole state, and has a small overlap 
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with the two low-lying particle configurations. 

There is no reason why the addition of a neutron to °%Cu should 
not affect the coupling between the proton and the °?Ni core. While 
this model assumes that the presence of the extra nucleons does not 
significantly perturb the core (and hence we can use the excitation 
energy of the first 2" core state for the phonon energy fw), this is 
only approximately true, since the addition of two neutrons changes 
the phonon energy from 1.172 MeV in °*Ni to 1.348 MeV in °*Ni. bie 
therefore seems reasonable to allow the proton coupling strength ae 
to vary somewhat from the value which gives the best fit to the levels 
of °*Cu. However, Heyde and Brussaard (He 68) found in their 
calculations for ‘*°La that the energy spectrum was relatively 
insensitive to the individual values of oD and eu? provided that the 
sum 3 - cn was a constant. It is doubtful therefore that allowing 
20 to vary would produce a much better fit. 

It is also considered unlikely that including N = 3 phonon states 
in the calculation would significantly alter the low-energy spectrum 
ore rca. Table IX shows that the lowest energy levels of ’*Cu are 
built mainly on the lowest six states of ®3cu, which are composed 
predominantly of O and 1 phonon states. The contributions from 
two phonon states are in many cases not negligible, but those from 


three phonon states almost certainly will be. 
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CHAPTER 6 


CONCLUSIONS 


The (p,ny) reaction has been shown to be an extremely powerful 
tool for studying nuclei in the nuclear p-f shell. Because of the 
fact that the reaction mechanism is well described by the Statistical 
Compound Nuclear model, angular distribution measurements alone are 
often sufficient to uniquely determine the spin of the decaying level 
in the residual nucleus, as well as the multipole mixing ratio of the 
Y ray. 

In certain instances, the angular distribution measurement does 
determine the level spin uniquely, but cannot distinguish between two 
possible values of the mixing ratio. Linear polarization measurements 
can often be helpful in resolving this ambiguity. Such measurements 
also provide information on the parity of the radiation, and hence on 
the relative parities of the initial and final nuclear states in the 
decay. 

When the angular distribution of the decay is nearly isotropic, 
little information can be deduced concerning the spins and mixing 
ratios. For the (p,ny) reaction near threshold, this situation is 
expected to occur whenever the spin of the decaying level is less than 
is the spin of the target nucleus. In 


or equal to | 1, where | 


Te T 


such cases, the measurement of a y-ray triple correlation can often 
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eliminate one or more of the possible spin assignments. 

In the present experiment, we have used all of the above techniques 
to obtain unique spin assignments for almost all of the levels in ACG 
below |] MeV excitation, and multipole mixing ratios for most of the 
Y rays. The spin assignments deduced in the present work are generally 
in good agreement with other available results. Values for the 
multipole mixing ratios determined in this experiment are also in 
fairly good agreement with previous works, although our values tend to 
show somewhat higher E2 contributions for many of the decays. 

From a theoretical standpoint, the Intermediate Coupling model 
given a fairly good pee cuesnice of many of the low-lying energy levels 
of ®*Cu, but there are still some aspects of this nucleus which are 
not clearly understood. The failure to predict the second _ state 
at 278 keV excitation must be regarded as one of the most serious 
failings of the model. It is possible that including mixing between 
particle and hole states in ®*Cu would yield a better overall agreement 
with the experimental energy spectrum. Calculation of properties 
which are more sensitive to the details of the nuclear wavefunctions, 
such as electromagnetic decay properties, would undoubtedly give 


further insight into the validity of this model for S Buy 
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[. INTRODUCTION 


The study of the electromagnetic decay of an excited nucleus is a 
widely used technique for the determination of the properties of nuclear 
states. The radiation pattern (i.e., angular distribution or double 
correlation) of the Y rays emitted by these states is intimately related 
to the multipole character of the radiation, and its measurement can 
provide useful information on the angular momenta of the nuclear states 
involved in the decay. Measurement of the double correlation can also 
provide information on the values of the multipole mixing ratios, which 
are proportional to the relative amplitudes of the different multipoles 
contributing to the decay. 

The information obtained from the measurement of a mixing ratio is, 
in principle, twofold, since both a sign and an absolute magnitude are 
required to fit the experimental data. Until quite recently, the sign 
of the mixing ratio has been of little importance in the comparison of 
experimental data with the predictions of nuclear models, since it was 
impossible to relate the phase conventions of different authors back to 
a common starting point. However, with the recent development of a 
phase consistent theory for y decay by Rose and Brink (1967) (hereafter 
referred to as RB), this problem has been largely overcome, and the 
measured sign of the mixing ratio has now become a useful parameter to 
be compared with model predictions. 

The amount of information obtainable from the analysis of a double 


correlation measurement is severely limited by the fact that the number 


105 


qestoun to esiitsqeig sdt Fo ncitenimessbh ors 407 supinrioed cae 
alduob to noitudiateib qelugts , 2.1) m19396q soiseibery sat .osede 


betelet yisismisni 2i 2atet2 Saens yd bassime eys7 Yond to (noiteh: 17109 


zetede iseloun end to sinsmom aslugns edz no noizamiotni fytseu eblvarq : 
cate nép noldetsy103.sliducb ad To snenieruassM .ysoeb s1t7 ni bevioval 
Adidw .20ise7 ontxim slogitfum eft to eaylev si3 no aoldemyotal eblveiq 
estogit tum dnsisttib of? to 2obud il amie svideley siz of Isnoisvoqotg 816 
.y6oeb orld of gntgudin3909 
21 Gigs) batxim s to Jnameivesem oz mov? benissdo noisemotnl eit 
sis sbhutinoem siuloeds né bas agile » dtod soniz ,blotows ,siqlonizq ab | 
fele sit ~yisnsos, stiup [ignl .s3sb fsinemiisqxs arid 21% of boviuper 
Yo foziwsqmo> di ni sofesiogml sistil io nasd asd oles paleie orl io "i 


26w 3] sonle ,2isbom veeloun! to 2n0igaibssg scx diviw 6seb isanamiteqxs 


07 46d eroflive SAeisTT 1b 3 enolinavnos sesiq siz e76fe1 03 aldlezoqmt " 
6 to signe a jng29 ent rtiw  vevewor sIniiog enlsisze nomnos § 


‘sertesreH) (rae!) Anii8 bas s20A yd ysasb y yo? yiosd? snez2ianos wae r 


yi all (SmO3718V0 ent nesd' zedimaitotg 2ina , (88 26 of bevieter 7 
a 


106 


of unknowns in the problem often exceeds the number of experimentally 
measurable parameters. This restriction can be overcome to a certain 
extent by measuring the circular and linear polarization distributions 
of the y decay. Since the linear polarization double correlation 
depends on the parity of the radiation (the polarization insensitive 
double correlation does not), one additional parameter is introduced. 

Further information about the decay can be obtained from the 
measurement of the angular triple correlation (two radiations in 
coincidence). The number of experimentally measurable quantities is 
greatly increased in this case, and (although more unknown parameters 
are again introduced) the problem is, in general, overdetermined. 
Linear polarization triple correlations can also provide information 
on the parities of the radiations. 

Since the measurement of triple correlations represents such a 
powerful tool in nuclear spectroscopy, it is clearly desirable to 
extend the phase consistent theory of RB to cover the more general 
case. The formula, which is derived in detail in RB, for the transi- 
tion amplitude for the y decay of an initial state | JaMi> to a final 
state | JoMo> is used as a starting point in this derivation. From 
this expression, formulae for the double correlation and polarization 
double correlations of a single Y decay will be derived. Only cases 
where the initial nuclear state possesses axial symmetry with respect 
to a space fixed z-axis will be considered, and the initial state Ji 
can then be completely specified by the relative populations of its 


magnetic substates (the population parameters P(M:)). 
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Time-dependent perturbation theory allows one to relate this 
amplitude to the transition amplitude for the emission of two 
successive Y rays, and the angular triple correlation and polar- 
ization triple correlation formulae are then derived in a phase- 
consistent manner. The formalism can be extended to include an 
arbitrary number of successive Y decays, and a recipe is provided 


for writing down the correlation in such cases. 


11. Y-RAY DOUBLE CORRELATIONS 


Consider the transition from a nuclear state |JiMi> toa 
state | JoM2> by the emission of a Y ray in the direction k 
with a given polarization é, Since the electromagnetic field is 
a transverse field, it is clear that (at most) two linearly 
independent vectors orthogonal to k are required to describe the 
polarization (which is conventionally taken to be parallel to the 
electric vector of the radiation). 

lt is advantageous to choose as a basis set, unit vectors 


which have convenient properties under rotations. A suitable set 


of basis vectors for describing the polarization is 


nw 


where e_ has been chosen parallel to the direction of propagation 


z 
> 
k. The vector with q=+ 1 represents a state of right circular 
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polarization, and that with q=- 1 a state of left circular polar- 
ization. 
The amplitude for the emission of aY ray with circular polariz- 


n~ n~ 


ation € = ea is given in RB 


E 


Mi -Moq (R) (2) 


1 
AMA Mp iii (se) 2 ¢ <M LT ay oe | JoMe> D 
where R= (9,0,0) is the rotation which takes the z-axis from the 
initial (space-fixed) direction to the direction of propagation k 
of the Y ray (see Fig.37). The definition of the rotation matrices 
used here is that of RB. 
The Tee Of (2) are the phase-defined electromagnetic transition 
operators of order L (explicit definitions are given in RB). They 


are defined such that they transform under time reversal according to 
aM ees Tr. 
Sweats 1 (3) 


Where 0 is the time reversal operator. This, together with the 


requirement that the nuclear wavefunctions transform according to 
ejum> = (-)9 | sy (4) 


e e e e <1> 
[5 suffictent to ensure that the matrix elements <JiMi|T)y | JoMo> 


are real. 
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The index 7 specifies the electric or magnetic nature of the 
radiation (wm = 0 for electric multipoles and 7 = 1 for magnetic 

A : AA oh . es oti 
multipoles). The parity of the radiation is then (-) ; 

In general, the orientation (i.e. the projection M2) of the 
final nuclear state is not observed with the result that the transi- 
tion probability is summed over all Mz. Moreover, the initial state 


generally exists in several magnetic substates M; with relative 


populations P(M,). The total transition probability is then given by 
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Since the matrix elements of Tis are real, the complex conjuga- 
tion on the matrix element can be dropped. 
The Wigner-Eckart theorem separates the dependence of the matrix 


elements on magnetic quantum numbers 
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and the two rotation matrices can be combined. 
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Mi-Mo- 
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Since nie (00620) a= P_ (cos®) , Eq. (5) becomes 
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Due to the assumed axial symmetry of the initial state |JiMi>, 
the transition probability is independent of the azimuthal angle 9. 


The sum over My in Eq. (8) can be evaluated explicitly to give 


(.)4 +Ji-Jot = Lt =k JinMy 


JIW(JaJiLL'3kd2) (-) (Jada3Mi-Mi{/kO0) (9) 
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which contain all of the information about the initial state in the 


y decay. Defining (RB) 
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In the angular distribution and polarization distribution experi- 
ments, the relative yield only is determined as a function of angle, 
and it is clear that only the ratios of the reduced matrix elements 

<T1> 


can be determined. For this purpose, multipole mixing ratios, Oy 5 


are defined 


te) th yA 
5 L va 


da [175 []422/ b 3) 


where (L 7) is the lowest allowed multipole in the decay J, > Jo. 
Since the reduced matrix elements are real, the mixing ratios are 
also real numbers. Normalizing Eq. (12) by dividing by the factor 

2 2 (ye aes and neglecting constants, the circular polarization 


double correlation is given by 
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If the circular polarization of the y ray is not measured, 
Eq. (14) must be summed incoherently over the circular polarization 
quantum number q, to arrive at the formula for the polarization- 
insensitive double correlation. 
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In deriving (16), the only restriction placed on the initial 
and final nuclear states was that they have definite angular 
momentum J, and J» respectively and that the initial state be 
formed in axial symmetry. If these states also have definite 
parity (and parity is conserved in the y decay) then Layee oe 


= 1. The sum over k _ is then restricted to even values only. The 


normalization of (16) is such that ap = 1. 
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The B, (Ja) coefficients contain all of the information about the 
initial state of the system, through the population parameters P(M1). 
There are three classes of population parameters which are often 


encountered in practical applications. 


Isotropy. P(Mi) = 1/512 for all Mi. Since the populations are 
equal only Bo(Ji) is non-zero, and the angular distribution is iso- 
tropic. This is to be expected, since the isotropy of the substate 
populations has effectively removed the axis of symmetry from the 


bineLibaies tate. 


Alignment. P(Mi) = P(-Mi). In this case, 


Be(J1) = 8a, 09 J, (-)7? (Jid1300|k0) P(0) (17) 
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Polarization. P(M,) # P(-M,). The By (Ji) coefficients are, 


in general, non-zero for both even and odd values of k. 


IV. CIRCULAR POLARIZATION DOUBLE CORRELATION 


Equation (14) gives the distribution of y rays with a given 


circular polarization, specified by the quantum number q. Following 
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the usual definition for the degree of circular polarization, 
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: nie ene wit (a) 
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are 


PEON atop de a SINT gu tg (18) 


This definition limits B (8) to lie between +1 and -1, with 
P_ (6) = 0 corresponding to an unpolarized y ray. 

lf a parity conserving transition occurs between two states of 
definite parity, the sum over k is restricted to odd values only. 
Unless the initial state is polarized in some way, as in experiments 
using polarized targets or polarized beams, the By (J2) vanish for 


odd k and there is no measurable circular polarization. 


V. LINEAR POLARIZATION DOUBLE CORRELATION 


Since the two unit vectors in Eq. (1) form a complete set of 
polarization states for the transverse electromagnetic field, a 
State of arbitrary linear polarization can be expressed as a linear 
superposition of the two circular polarization states, and the 
transition amplitude for the emission of a linearly polarized photon 
can be expressed as the same superposition of the two circular 
polarization amplitudes. Specifically, the transition amplitude for 


the emission of a photon linearly polarized along the xdirection is 
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Notice that the coordinate system is chosen such that the 
original z-axis and the direction of the y ray (which are connected 
by the rotation R = (¢,8,0)) define the x-z plane, and Eq. (19) 
then gives the amplitude for the emission of a photon linearly 
polarized parallel to the reaction plane. The amplitude for the 
emission of a Y ray linearly polarized at angle € to the reaction 
plane is found by rotating the coordinate system through angle & 
about the final y ray direction, and is then given by Eq. (19) 


with R= (¢,6,&). The transition probability is 
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and the linear polarization double correlation is given by 
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By analogy with the case for circular polarization, the degree 


of linear polarization is defined as 
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Note that parity considerations again restrict the sum over k_ to 


even values. 


VI. Y-RAY TRIPLE CORRELATIONS 


Consider the process shown in Fig.38, where two y rays are 
emitted sequentially. Using first-order time dependent perturba- 
tion theory, it can be shown (Appendix 1) that the amplitude for 
this two step process is just the product of the amplitudes for 
the individual decays. In general, the cascade decay |J1M1> > 
| oM2> S2 | 3M3> can proceed through any one of the (2J2 + 1) 
degenerate states |JoMo> » and the total amplitude is given by 


a sum over the intermediate states 
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since k,; + ko + k3 is an even number (Appendix 2). 
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This is the basic formula for all triple correlations. 
[lf the circular polarization is not measured for either of the 


Y rays, an incoherent sum over gq; and qo gives the formula for 
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the direction-direction triple correlation 
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Parity conservation in the y decay, plus the requirement that all 
nuclear states have definite parity, limits the sums over kp and k3 


(and therefore k,; also) to even values. 
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Eq. (31) over all angles {2 and sum over the circular polarization 
quantum number q2. Since this process is then equivalent to a 
single transition | JiM> ma | JoMo> the general formula should reduce 
to Eq. (16). 
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After the sum over qo (= +1), this reduces to Eq. (16) for 
a y-ray direction double correlation, apart from a normalization 
PacCtOneOt medi. 
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Since the UL (L3Ji3J2) are less than unity in absolute 
value for k,; #0, the effect of the unobserved transition is 


to decrease the anisotropy of the double correlation. 


VItt. TRIPLE CORRELATIONS INVOLVING CIRCULAR POLARIZATION 


[lf the circular polarization of y, is measured and that of yz 
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The circular-direction triple correlation is then given by 
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where the superscript (1) indicates that the polarization of Yi is 
measured. Parity restrictions require that k» be odd. 


If it is the circular polarization of y. that is observed, a 


sum of Eq. (31) over qi obtains the formula for the direction-circular 


triple correlation. The direction-circular triple correlation in this 


case is 


23 


’ , 
- Tal : 
a es ov LP 
Cid ey Tey us 4 (rt) 8 t= ( a, hee 
ied casa ae 
"ig 
= tot! cbeste dt 7 ; ; 
‘ pet git’ s s £ (<) +1} (Nsdecd) A q x 
(i) te! 


; ~ a 
(se 19) yo (elsb'etst) 7 (gbeb' 1dr) credit : 


: oe 
Vd navio Aedi ef noltefatios siqing noitas1ib-r6fvatt Msi j 


i ui 


= 


( oe, See CRETE ane iad 


aie ented 1m 
. ‘wd 


- : 
peer h smth itera ye Ty (ete cdyA f Ucbrb tay gy 


gtd 
Vor sd 
- . 
¢ F : i 
(se) (<2. 4) aled oh (ebst sdsd) , os 


ong 


. : Fag kt oe 
2i «y¥ to noltesixelog edt ferns eets60ibni (tr) Iqi We TOUE MEI e aft 


wbbo ed gt Asi sviups. enolaaiyseen wine 
‘s sbaviszdo- ei ol id to nolaaxt+6leg asivotts > at 2 ale 
‘stumed te enintdo 1. revert) a -pa 


~ 


ao 
Site 


Pray © 


qelpo7, ioraol 43971 ei y 410 
pie nb aoiseloroo ¢ gin: bs ystlabbandatn 

“7 4" 7 
a firm ae 


- - 
4 


p (2) 


(21,22) = bor Bed Jadee ) Yohe 2 fbag by) 
Cc W(21 , 22) a ka ie aera 


Lana! 


{1+ cay brtMi tls 4a Koy, 


i] 
kikoks (bibs JiJd2) 


“ r% Ae ay Rea rs ena 
Loto 
Lo Mo! 


i] 
x RL, (5 lo ded) Pet (ae Glo) (43) 


For the case where the circular polarization of both y rays is 


measured, the circular-circular triple correlation is 
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[X. CORRELATIONS INVOLVING LINEAR POLARIZATION 


To derive the triple correlation where the linear polarization of 
one (or both) of the y rays is measured, the corresponding transition 


amplitude (s) must be expressed as a linear superposition of the two 
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circular polarization states. For the case where the first Y ray is 


linearly polarized at angle — to the reaction plane, the transition 


amplitude is given by 


a fa S > 
Anny Erk) = 2 Aye, (Ei) Ay th (he) (45) 
2 
1 
(kik)? tee art 1 <> 
qd, JiMi|T uy. |J2Meo> 
2m ¥2— Mo qi Lit, Lote PUT a 
Ly Lo R ) 
x Mi -Mo qy (Ri) q2 es) 2Mo|T 0” -y, [YaMs> Dy? -M sq, ( 2 
with Ry = ($1,61,81) and Rp = ($2,02,0). 


Proceeding in the usual way 


gestions +1 1, '+1 
Ce A(Ly by") ean ape 


Wed (2, So), = y ah 7} 


ky Lim qiqi' 
Peak 


1 

Same aca can (Costs! ae ray ye a eneneUi 2) 
Loto Koike 
Lo't2' 


RET? (lolotdoda)vP, Ud OtOp iy EieO> 00a) (he) 


Nee | ea he ke 


where 


etn (Lib 'JiJ2) (-) 41h eco d ieee (47) 


Jia hy Jo 
Ke yee dr-dqe Kodi qi iy el 
ka) Ko k3 


i 
a et ¥ tent? sds. wet! Rae add 109 oe aoitest 


notatensys atit ea alls nolioset adi os 3 steno je si 
yd nsvip: ai abut 


p - 
(2#) (3) SPA COD yin i = (aa) 
I | 
- m, . Se siu 
cit hig hatha kk 2 
git- ki iP eWel att 4p se “S\ ars 
of) oJ ot? gst. a tJ } : 
'e we deel i _JPatst ak CD 5 oe 1H 
(0,06, +6) = oh tas (63.20.10) ;" * 


yew |sueu eat ni gnibess 


i+? pm. feat ‘rte td) oO 7 Da 
cai ad | Citints (Cet) ae i = (ie DP 
i c) 
Ts My { 


(obab)sd) ea {0 ctda $ 
eng ets ry 


ja ; o o'rp7 ep SP 
(o#) (0,62. 23,19) etod t (ebeb' sted) ont x 


Al ga 2 4d Mid ry Aebabtededd | i 
— > fe ,? 7 


pares Ard bari 7 p an ; 
- : 


7 
7 


and 


p 41-q1' qo-q2! 


kikok3 (81,22) 


lf the direction only of Yo 


We (1,22) = } By 1) ) 


lati 


126 


= kok } (kok33K-K|k10) 
ioe 
Nise Meee GH eer, Ur CO ACEy 
is observed, a sum over q2 gives 
Tit) 7, '+1 


aide es 
2 ; 1 
91qi 


Ly ara! 


so) 
Loto 
Lo ‘'t,' 


Ato) 


» {4 “ (2) pene cae alas 
(49) 


koks 


; Lee 
eas (Lita Jad2) Ry, tLab2'Jods) pated? (omen ss0) 


Kikeks 


kikok3 


and the linear polarization-direction triple correlation is given by 
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Lf it is Y2 which is linearly polarized at angle & to the 
reaction plane, the transition amplitude is a linear superposition 
of the amplitudes for the two circular polarization states qo 


Proceeding as above, 
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As before, a sum over q, is required if the polarization of y, 


is not measured. The direction-linear-polarization correlation is 
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X. AN ARBITRARY NUMBER OF TRANSITIONS 


The extension of the formalism for an arbitrary number of y-ray 
transitions can be inferred from the study of three successive y- 


transitions. The amplitude for this process is given by 
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The quadruple correlation, for the case where the circular polarization 


of all three y rays is observed, is given by 
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By comparing this expression with Eq. (14) for a single y-ray and 
Eq. (31) for two Y rays, a program for writing down the correlation 


for an arbitrary number of transitions can be detailed. 


(1) Beginning with the initial state | a> in the sequence, write 
a different k. value for every nuclear level and every y ray 
in the sequence, except for the last y ray and the final nuclear 
state. Also if any of the Y transitions are unobserved, do not 
write a k. value for this transition, and make the k. values 
for the two nuclear states connected by this unobserved transi- 
tion the same. 

(2) tn the correlation formula, write down and sum over all of these 
k. values, and insert the By Wa) coefficient for the initial 
State. 

(3) For every observed y transition, write down the factor A(L,,L.') 
and sum over L.; Tes L.', TT. 

(4) For every y decay (except the last one) which connects the state 
[J «> to the state ay by multipoles L., L.', we insert a 


factor which depends on what we observe about the y ray. 
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iii) {f the direction only is observed, we insert the factor 
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iv) if the y-ray is unobserved, we insert the factor 


For the last observed y-ray in the sequence, the above rules hold 
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The angular functions are defined as follows: 
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proceeding towards the first, insert the following factor for 
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For the first two D-matrices in this sequence, write a Clebsch- 
Gordan coefficient (kk';«K'|k''K"') to couple k for the 
second last y decay to k' for the last decay (with the 
corresponding Kk from the D-matrices) to k'' for the initial 
state in the second last decay. For each remaining D-matrix, 


write a Clebsch-Gordan coefficient to couple the k of the 


D-matrix to the last k in the preceding Clebsch-Gordan, to 
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give the k for the initial state of the transition. For the 
last Clebsch in the series, set the projection of the final k 
equal to zero. 

(3) For each Clebsch (kk's« K'|k''K+K') insert the factor 


k k' | k'" and finally sum over all the kK. 


After completing this procedure, we are left with the angular 
correlation formula W(Q21,22, . . .Q,). If polarizations are being 
measured, the polarization correlations P(Q, . . . %,) can be 


defined in the usual way. 
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Figure 37. The Euler angles (a,8,y) of the rotation R which 
takes the coordinate system S(x,y,z) to S'(x',y',z'). 
This transformation is accomplished by a rotation 
of angle a about the z axis (to S''), followed by 
a rotation of angle B about the y'! axis (to $'!'), 
followed by a rotation of angle y about the 2z''' 
axis (to S'). 
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Figure 38. A schematic diagram of the y ray cascade 
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APPENDIX | Y-RAY ANGULAR CORRELATIONS 


Consider the situation shown schematically in fig. 2, in which 
two Y rays are emitted sequentially. Using time-dependent perturba- 
tion theory, we can relate the amplitude for this process to that 
for a single y decay. 

We assume that the three nuclear states v; = [J M.> Gi=len zens) 
are eigenstates of a time-independent Hamiltonian Ko 


Ht, My ee (Al-1) 


The total Hamiltonian for this system is H = Ho +fiyy (t) where 
A yr (t) is a time dependent interaction which induces transitions 
between the eigenstates - The total wave function ¥ is a solu- 


tion of the Schrodinger equation 
Hy = any (Al-2) 


We assume further that the states vy; form a complete ortho- 
normal set for the system, and we can therefore expand the total 


wave function as 
Y(t) = ) alee (Ai-3) 
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The square of the expansion coefficients a, (t) gives the probability 
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of finding the system inthe state p, at time t. 
Substituting Eq. (3) into Eq. (2), and making use of the 
orthonormality of the wavefunctions vy; when integrated over. all 


space coordinates, we arrive at the set of coupled differential 


i ,(t) = ar(efly, + ag(e) Mie ef 2° + agle) fy, of rtwa)t 

Toes ahelmiee 2. 471i) pndcmn Game REN eeahe) ee 2 ee (ar-t 
testes Beet uta ee) see) Pens 25 estes 
where 

ne awe: 

Wo = (Eo-E3)/h 

Aig = ei awe 1 ¥5> a 


The terms Involving His sich describe transitions directly 
from the state y~, to the state 3, and while such transitions may 
in fact occur, we are interested only in the cascade decay, and 
these terms may be excluded. 

In order to solve the set of Eqs. (4), we assume that the 
interaction is turned on at t = 0, at which time the system is 


entirely in the state ~,;. We then have the boundary conditions 
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We also assume that for a short time after t = 0, we can neglect the 
time dependence of the a. (t) on the right hand side of these equations. 


The second of these equations is then 
ih a (t) = H., (t) aoe (Al-7) 


heene is a slowly varying function of time, we replace it with 


its value at t = 0, and can integrate (7) to get 
a(t) = ato) portmt _ 43 (al-8) 
With 
Substituting this expression into the third of the Eq. (4) gives 
ih a, (t) = fen (o) Hap (t) a2 femele - 1] (A1-9) 
1 


which can be solved in a similar manner to give 


-i(wytwo)t _ ipl 
a3(t) = (0) DR lie se a le Se at (A1-10) 
W, #2 Wy, + Wo Wo 
Since we have explicitly removed any contribution from the 
direct transition ~, + ~3, the population of the state 3 at any 


time t is a measure of the probability that the cascade transition 


has taken place. We therefore consider the term 
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Finally, we integrate (11) over all energies W, and Wo. of the 
photon states. This converts the second factor in the above to 


27™t, and the transition rate is 
d 2 
OS Goer eee 
2 
5 Ao (0) H 5 (0) | (Al-12) 
apart from constant factors. The basic result is that the amplitude 


for the two-step process is just the product of the amplitudes for 


the individual decays. 
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APPENDIX 2 k SELECTION RULES 


Consider the term in the triple correlation formula 


wie Sat Reig, Canc rey, (A2-1) 
Since Lj and L,' span the same range of values, this term cannot 

change when we interchange L; and L;'. Making use of the well- 

known symmetry relations for the 9-j symbols and Clebsch-Gordan 
coefficients appearing in Eq. (30), we see that such an interchange 
kotk3-ky 


introduces only a phase (-) , and we can then demand that 


k; + kp + k3 be an even number. 
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A DESCRIPTION OF CALCULATIONS IN THE 


INTERMEDIATE COUPLING MODEL 
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TiTLE= | CUPPLER) 
TPE eMail 
SOURCE LANGUAGE: FORTRAN IV (6G) 
ORIGIN: P.W. Green 


Nuclear Research Centre 
The University of Alberta 


ABSTRACT: The purpose of this computer program is to calculate the 
energy eigenvalues and eigenfunctions, using the framework of 
the Intermediate Coupling model, for either a single particle 
coupled to the surface vibrations of an even-even core (for 
odd-A nuclei), or for two particles coupled to such a core 


(for even-even or odd-odd nuclei). 


| ONE-PARTICLE CALCULATIONS - THEORETICAL FORMULATION 
For a system composed of a single particle coupled to a collective 


core, the total Hamiltonian of the system can be expressed in the form 


H. =H. +H +H (1) 


i G p INT 


where He is the Hamiltonian for the core, ts is the particle 


Hamiltonian, and H represents the interaction energy between the 


INT 


two systems. 
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The core is assumed to be capable of performing quadrupole 
oscillations, such that the nuclear surface can be described by 
R(8,$) = Roll + } a* Yo" (6,4)] (2) 
u 
where the afte the dynamical variables of the collective motion. 
Provided that the deviations from a spherical shape are small, the 


collective Hamiltonian can be expressed in the form (Ch 54) 


He =) fBla,|? + clo,|?} (3) 


u 


which is the Hamiltonian of a system of harmonic oscillators with 


frequency 
w= ()* (4) 


The parameters B and C in (3) are, respectively, the collective mass 
and the nuclear deformability. 
lt is convenient to rewrite the core Hamiltonian (3) in terms of 


creation and annihilation operators (b t 


y and b) for spin-2 phonons, 


which are defined by the equations 


I 


“i 
b= {Bun + in J 
i ts notliert oll 
a a Bae [Bus , lis ] (5) 


° si 
here 1 = ‘BO, -~. 
where u i 


The core Hamiltonian then becomes 


it 1 
fhe } (b bi + >) tw (7) 
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and has, as its eigenvectors, a series of states with energies 
Ey = (N+ 2) tw (8) 
N 2 


N can be interpreted as the number of phonons associated with the 
particular state. 

For a given single particle orbital, the particle Hamiltonian 
th is a constant, having as its eigenvalues the single particle 
energies of the average spherical potential. 

The interaction Hamiltonian represents the energy due to the 
coupling between the extra particle and the quadrupole oscillations 


of the core. To first order in the Sis it can be written 


die ae eel Yo" (6,4) (9) 


u 
where (r,9,6) are the co-ordinates of the particle. 
In calculating the matrix elements of (9), we will need the 


value of the radial matrix element 
k = <n'2'[k(r) [n> (10) 


which will depend, in general, on the quantum numbers n and Q describing 

the initial and final states of the particle. The dependence of k 

on these quantum numbers is not critical, however, and it will be 

assumed that k is a constant. Again, it is convenient to write the 

Hamiltonian in terms of phonon creation and annihilation operators, 
Hiyt = 7 CD)? toe ETB aE)e bs Tye (ee) i) 
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dimensionless coupling parameter 


z oe 
We shall work in an angular momentum representation. The basis 


States of the system are chosen as those of the non-interacting system 
[NRj;tM> = ) (Rj;M-m m| IM) [NR M-m> | jn> (13) 
m 


representing a particle of angular momentum j and z-projection m, 
coupled to a core state characterized by N phonons coupled to a core 
spin R, to give a state of total angular momentum | and projection M. 

The eigenfunctions of the total Hamiltonian (1) will be linear 
combinations of these basis states, of the form 

a 
[E,;Im = } A (NRJ§L) | NRJ; 1M> (14) 

The expansion coefficients A (NRj!) and the corresponding energy 
eigenvalues in Eq. (14) can be found by calculating the matrix elements 
of Hos and diagonalizing the resulting matrix. In this representation, 
both the core Hamiltonian H. and the particle Hamiltonian Ms have only 


C 


diagonal matrix elements 


<NRJ;IM[H, + H[NRJsIM> = (N + 2 tw te, (15) 


where hw is the phonon energy and ej the single particle energy for 


the particle in orbital j. 
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The matrix elements of the interaction Hamiltonian are 


Les PB hea ea tie ai Shig pelle 


= - (Bhwe 2 (Raj sgMtem! m | 1M) (RjsMem mt) (=)? 
mm' 


x <NIRIM' =n! [bt eae b | NRM=m> <j'm'|¥> |[jm> (16) 


The Wigner-Eckart theorem* can be used to separate the dependence 


of the matrix elements on the magnetic quantum numbers 

<NIR'M' =m! [b_T [NR M-m> = & (R23;M-m -1[R'M'-m') 

R! 

x <n'R'||b' | |NR> 
<j! 'iyb (6 a fe asl (j A lj? ) <j' lly iitee 
jim! |¥2 (8,) |jm> = mw (j2;mujjim') <j 21 Ij 

ji! 

<N'R'M' =m! |b) [NR M-m> = <NR Nm|b, [N'R! M'=m!'> (17) 


“nw 


where j = ej shes The last line in the above is required because the 
Wigner-Eckart theorem can be used only with irreducible tensor 
operators. The creation operators b are irreducible tensors, but 


the annihilation operators b. are not. 


The definition of the Wigner-Eckart theorem used here is that of 
de Shalit andsiaimi. (Sh 63), 


<J'M! |T S| uM> aia! (Jk3Mc|J'M') <y'[|T, | [ue 
k J! k 


Because of this definition, our reduced matrix elements are a 


factor of J' larger than those of Choudhury (Ch 54) and Choudhury 
and O'Dwyer (Ch 67). 
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Using the techniques of Racah algebra, the sum of four Clebsch- 


Gordan coefficients reduces to 


Deen ren VSAM (ROMER Ie 12) 0h 3) nemee agnG. 3 


and the complete term for H is 


INT 


= Be mag YIN je yy| [se tne" | [BE | |e 


(pick <NR| |b" | [N'R!>} W(R'f'Rj312)6 6 


4 iy) Omm! (19) 


It is now necessary to calculate the reduced matrix elements 
appearing in (19). The results for the spherical harmonics are given 


by several authors. We use that of de Shalit and Talmi (Sh 63) 


L+Q! 
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To calculate the reduced matrix elements of the phonon creation 
Pro cee : ‘ sn 
operators bY » It iS convenient to introduce the coefficients of 
fractional parentage (cfp's) for an N-particle system. Using these 


coefficients, an N-phonon state can be expanded 
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where [2 (ai7R1) 2 R|} 2 aR] is the cfp, Ri and R are the total 


angular momenta of the (N-1) and N particle systems respectively, and 


&, and a represent additional quantum numbers needed to completely 
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specify these states. We are concerned only with states for N < 3, 

and since these states are completely specified by the number of phonons 
N and the total angular momentum R, the quantum numbers a and a, will 

be suppressed. 


Acting on this state with the annihilation operator by gives 


b, [NRM> ie) (oN tie.) 2R | } 2NR] Y (Ry2;m M-m|RM) 
Ry m 


x |N-1 Ry m> Sen 


wT) (24! (a1) 2R| 32 \R] (R 2;M-u RM) |[N-1 Ry M-p> (22) 
Ri 


where the factor YN comes from the fact that, for the annihilation 


operator acting on an N-particle state (Dw 67) 
b,  |N> = YN |N-1> (23) 
We then consider the matrix element 


<NIRIM! |b, |NRM> = <vRM |b, T|N'R'M'> 


= W JY (27! (Ry) 2R]} 2"R] (Ri2s-u ulRM) 
Ri 
Diep Sl Tre Sh ye 
Poe) toe N . 
= YN[2 (R') 2R|}2 R] (R'2;M-p HRM) 5 (24) 
N',N-1-M!'M-y 
Using the Wigner-Eckart theorem gives 
<NRM |b, T|N'R'M'> = 4 (R'25M"u | RM) <NR| |b" | [N'R'> (25) 
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and comparing (24) and (25), 
<NR| |b’ |N'R'> = /N Rl 2hwe (R') 2R|} 2NR] (26) 
A list of the reduced matrix elements for N < 3 is given in Table XII. 


1!i TWO PARTICLE CALCULATIONS - THEORETICAL FORMULATION 


For a system composed of two particles coupled to a vibrating core, 


the total Hamiltonian is 


H =H. + H (1) 4+°H* (2) + 4H + H 


~ C iy (2) + Hi2 (27) 


INT 


The notation is the same as that used in Sect. |, with the 
indices |] and 2 referring to the first and second particles 
respectively. The term Hi2 is the residual interaction between the 
two extra particles. 

It is possible to diagonalize the total Hamiltonian (27) directly, 
in a manner analogous to that given in Sect. |, to obtain the eigen- 
values and eigenvectors of the composite system. This procedure is 
often impractical, however, since it requires the diagonalization of 
very large matrices. 

The sizes of the matrices involved can be substantially reduced 
if the problem is divided into two steps. One first finds the 
eigenfunctions for one of the particles (say, particle 1) coupled to 


the core, by diagonalizing the Hamiltonian 


H, =H, + HM) + H (28) 
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lt is reasonable to assume that those eigenfunctions of H, which lie 
very high in energy will not substantially affect the low-energy 
properties of the core plus two particle system, so one may then 
choose a subset of these eigenfunctions as a basis set for the two 
particle calculation, and diagonalize the total Hamiltonian to obtain 
the results for the complete system. This is the approach which we 
follow in this work. 
The basis vectors for the two particle system are 
JordijostM> = ) (Jajo3M-m2 m2|tM) |oydi M-mg> |j2m2> (29) 
M2 
where a, iS an additional quantum number to specify which state of 
spin J, we are considering. The states Jordi M-m2>, which are obtained 
by diagonalizing the Hamiltonian (28), are of the form 
JayJy M-m,>= ) ee (NRJ adi) [NRJisJa M-m2> (30) 
NRjq 
The particle Hamiltonian H (2) is diagonal in this representation, and 


has as its matrix elements the single particle energies for particle 2. 
<oad aj 23 1M|H, (2) [ordaj 25 1M> = €2 (31) 


H, is also diagonal, and has matrix elements equal to the energy of 
the state lord; Ma Maer. The only off-diagonal terms come from 


Hi yy (2) and H,5. 


Considering H (2) + f.inSteawenhave 
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The matrix element containing the boson creation and annihilation 


operators is expanded in terms of its component states to give 
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x <NIR'M' =m, ‘=m, "| (-)"b_ + b [NR M-m,-m,> (33) 


What follows now is just Racah algebra. The Wigner-Eckart theorem 
is again used to separate the dependence of the matrix elements on 


magnetic quantum numbers, and the sum of six Clebsch-Gordan coefficients 


reduces to two Racah coefficients. The final result for H ayy (2) is 
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This term is diagonal in the co-ordinates of the first particle 
(j1), since the interaction does not involve particle 1. 
Finally, we consider the residual interaction between the 


particles Hio. 
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which is diagonal in the core states |NR M-m2-mi>. The interaction Hi2 
must be a scalar and it is reasonable to write it as a sum of scalar 


products of two irreducible tensor operators, of the form 


Hige= ) Sov, .(-) "Stes G1) at ee (2) (36) 
k k k 
kK 
where the operator T, “(1) acts only on the co-ordinates of particle 


k 


1, and ae(2) acts only on the co-ordinates of particle 2. The 
parameter Ve describes the strength of the interaction. The matrix 


element in (35) can then be written 
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As before, the sum of six Clebsch-Gordan coefficients can be 


reduced using Racah algebra, and the result for H,. is 
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x <j htt 2) PY 2? er) | W(kjaida'R3ja'di) 

MEW (dik I joldis 15) (38) 

A common choice for the residual interaction is a 6-function 
Hi2 = Vo 6(ri-ra) (39) 


It is more convenient to calculate the matrix elements of the 6-function 
interaction directly, since the multipole expansion contains an infinite 


number of terms. The result for Hio is 
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Fo = dr 45 Rani (r) R (r) R (r) R (r) (40) 


no! Xo! nih, n2Xo 


We again assume that the radial integral Fo is independent of the 
quantum numbers n and & and the strength of the 6-function interaction 
is described by the parameter 


N= VoFo (41) 


111 PROGRAM DESCRIPTION 


A) ONE PARTICLE CALCULATIONS 

For one particle calculations, the parameters of the model are the 
phonon energy (fw), the maximum number of phonons N, the single particle 
orbitals (2,j), the single particle energies (e5), the total spin (1) 


and the coupling strength (&). 
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For each set of input data, the program first sets up a matrix 
describing all possible basis states, composed of all core states with 
N phonons or fewer, coupled to one of the single particle orbitals. 
This matrix is printed, with an identification number | for each 
basis state. For each total spin value |, the program scans the list 
of all basis states, and picks out those which can give rise to a 
State of spin lI. Using this subset of basis states, the program then 
calculates the diagonal matrix elements, and the off-diagonal matrix 
elements (19) for € = 1.0. Finally, for each value of the inter- 
action strength €, the program sets up the total Hamiltonian matrix, 
and diagonalizes it to find the eigenvalues and eigenvectors. 

Output from the program begins with the matrix of basis states 
mentioned above. For each total spin value, a list of the 
identification numbers (1) of the subset of basis states used for the 
matrix is printed. Finally, for each § value, the eigenvalues are 
printed, along with the expansion coefficients A™ (NRj 1) Ofpeaes GLH). 
if desired. The format of the output is self-explanatory. 

The program also has the option of writing the results on unit 7, 
which is normally the card punch. This allows the results of a 
calculation to be used as input to another program. LG otnis option 
is in effect, every calculation (specified by a total spin | and a 
coupling strength &) is assigned an identification number KM, which 


is printed just before the eigenvalues. 


B) TWO PARTICLE CALCULATIONS 
For two particles, the parameters describing the calculation are 


the single particle orbitals and energies of particle 2, the coupling 
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strength ~,, the total spin 1, the form and parameters of the two 
particle residual interaction, and the energies and expansion coeffic- 
ients of the one-particle states (14), which in turn depend on the 
parameters listed in Sect. A. 

The parameters of the one-particle states can be obtained from the 
punched output (on unit 7) from the one-particle calculation. The 
format is such that no rearrangement of these cards is required, and 
they may be used exactly as they are for input. 

lt is generally not necessary to use all of the one-particle 
eigenstates as basis states. Those which are required are indicated 
in the following manner. 

One first indicates to the program how many data sets (each of 
which contains all of the eigenvectors for a given spin and &) it will 
be required to search. For each of these data sets, one then indicates 
the identification number (KM) for that data set, and how many (NBASE) 
states are to be included. The lowest (in energy) NBASE eigenvectors 
for this spin will be used. The only restrictions are that the list 
of identification numbers KM must be in increasing order (since it is 
not possible to backspace the input unit), and that the total number of 
one-particle states (ZNBASE) must not exceed 40. Note that the eigen- 
values and eigenvectors are read from input unit 4. 

After all the data has been input, the calculations proceed in the 
Same manner as for the one-particle case, and the output is identical 


in format. 
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1V PROGRAM INPUT 


Excluding the input from unit 4 discussed above, all other input 
to the program comes from unit 5 (generally the card reader). The 


input cards are in the form 


IC, (DUM(J), J=1,7) FORMAT (12,8X,7F10.0). 


iC is a (two digit) control integer, specifying the meaning to be 
attached to the dummy variables DUM, and perhaps also specifying the 
data to be entered on the succeeding cards. Table XIII contains a 
list of the options. 

All of the input data is remembered by the program. That is, all 
of the data associated with a certain value of IC, once it has been 
entered, will be used in all further calculations, until another card 
with the same value of IC is read. 

All input data is copied to the line printer. At the end of each 
input data set (indicated by an IC = 10 card), the current values of 
all relevant parameters are checked for obvious errors. lf any errors 
are found, the calculations are aborted, and the program returns to 
get more input. This process continues until no errors are found in 


the current data values, at which time calculations proceed. 


Verne le FUESHINGS 


1) Although it is possible to make unit 7 the card punch and unit 4 
the card reader, this procedure is not recommended. lt is easier 


(and safer) to assign these units to a magnetic tape or disc file. 
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2) For the two particle residual interaction (Eq.(36)), the form of 
the tensor operators TQ) and T, (2) has not been specified in the 
current version of the program. This can be done by modifying the 
function subprograms TIK(K) and T2K(K) such that for a given value of 
K, TIK returns the value of the reduced matrix element 


He HT MA) | Ge and T2K returns the value <ja'||T, (2) [|j2>- 
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TABLE XII 


Reduced matrix elements for the boson 
creation operators. 


eee ae 


+ 2 

N R N! R! (<NR|| b | |N'R'>) 
1 2 0 0 5 

2 0 1 2 2 

5 2 1 2 10 

2 F 1 2 18 

3 2 2 0 i 

3 0 2 9) 3 

3 2 2 2 20/7 

3 3 2 Z NS 

3 4 2 Z 99/7 

3 2 2 h 36/5 

3 3 2 h k 6 

3 i z 4 90/7 

3 6 2 hk 39 


ny 
7" 


The presence of an asterisk preceding a number indicates that the 
reduced matrix element is given by the negative square root of that 
number. 
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TABLE XI11 


tinout) Cards for GUPPLE = 71 


DUM (J) NAME 
J=? 
1 IWHICH 
2 NFONON 
3 EFONON 
Z NSET 
1 1D 
2 NBASE 


REMARKS 


Causes program to exit 
ignored 


The remainder of this card is ignored. 
The program then expects a card 
containing a descriptive TITLE in (20A4) 
format. 


IWHICH | for one particle calculations 
IWHICH = 2 for two particle calculations 
FOR _IWHICH = | 


Maximum number of phonons 
0 < NFONON < 3 


The phonon energy fw (MeV) 
FOR IWHICH = 2 


The number of data sets (each containing 
the eigenvalues and eigenfunctions 
for core states with a given Ji) which 
will be read from unit 4. 


The program then expects NSET data cards, 
each of which must have IC=0, with the 
following information on each card. 


The data set ID number (KM) 
The number of states of spin Ja, to be 
included as core states. The lowest 


energy NBASE states will] be used. 


Restrictions 
lnbteeN SEE S10 


Y NBASE. < 40 
i 
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03 


04 


05 


NAME 


NPART 


N1(1) 

PARTL (1) 
PARTJ (1) 
PARTE (1) 


XIMIN 


XISTP 


X 1 MAX 


LPUNCH 


CEUEL 


LCOEFF 


KMSTRT 


MAXDIA 


MAXRNK 
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The number of single particle orbitals 
allowed for the particle (1 < NPART < 4) 
The program then expects NPART data cards, 
each with IC=0, containing the following 
information for each orbital. 


The radial quantum number ns 


The orbital angular momentum L. 


I+ 
ur 


The total angular momentum J; = 2 
The single particle energy e, 


The minimum value for the coupling strength 


e 
The increment in & 
The maximum value of & 


Calculations for each spin will be done for 
values of € from XIMIN to XIMAX, in incre- 
ments of XISTP. For one value of € only, 
make XISTP = 0, and XIMIN = XIMAX 


=] Results will be punched on cards 
=0 Punched output will be suppressed 


Should be 0 always. LFULL = 1,2 provides 
two levels of full print, which are 
useful for checking out the program, 
but should not be used otherwise. 


=] Print (and punch if LPUNCH = 1) both 
the eigenvalues and expansion 
coefficients of the eigenstates. 

=0 Print (and punch) only eigenvalues 


Only needed if LPUNCH = 1. The starting ID 
number for punched output. Subsequent 
data sets are numbered sequentially. 


Maximum number of diagonalizations 


The maximum matrix rank (< 50) 


a ¥ a sone aay v 
A Hoes 407 nolwvar ss hee 


a sills alien. taoen a 


;2 mutnemom 1eC¥enB taatdto aAT (1) 3THAS 


2 + ,4 = 1 mutnaniom ssfupne tesos eT 

;2 Kprens sfaiiyeq atenta aAT 

dtprsize enilquos oni ot suléev mumirin oT 
é) 


2 of gnemsyoni onT 
2 to sulev mimixem oAT 


jot snob ad. f iw nige dase tc? anoissfuated 


-otoni nl. X% of WIMIX “3°10 2evtev 
ino M t 9 tev Saini joi .dfe2tk Fo 24pen 
% bne 10 ‘= W2ix sAem 


ons ro bodonua sd Filw esivesh i= 
bezgearqque od {fiw suqzu0 berlonus O- 


pis 272 eyeb 


aa Yo yedmur mun eh 
(62) inex Xi9 98m aiumixe OAT 


(1) 


(1) LTAAS 
(1) 3FaAS 


06 


07 


08 


09 


KLIMIT 


VPN (1) 


VPN (2) 


NSPIN 


RSPIN(1) 


RSPIN (2) 


XDELT 


REMARKS 


For IWHICH = 2 only. The two particle 
residual interaction is taken as a 
multipole expansion, and terms up to 
and including k = KLIM!IT are to be 
considered. OFS KEIM <a: 


V, for k = 0 


VGtrorak 
etc. 


] 


Number of total spin value |! for which 
calculations will be done. 


1 < NSPIN < 6 
First I-value 


Second I-value 
etc. 


For IWHICH = 2 only. The two particle 
interaction is taken as a delta 
function. This entry is the 
parameter n = VoFo describing the 
Strength of this interaction. 


Not used 
End of input. Check the current value 


of all parameters, and if no errors 
are found, begin calculation. 
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